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1.0 INTRODUCTION 


In Ref. 1 a Solar Power Satellite Transmission System incorporat- 
ing automatic beam forming, steering and phase control was presented. 

The phase control concept, which was discussed in great detail in 
Ref. 1, centered around the notation of an active retrod i recti ve 
phased-array , ([Refs. 2,3]. Figure 1.1 illustrates the major elements 
required in the operation of an SPS system which employs retrcdi recti vity 
as a means of automatically pointing the beam to the appropriate spot 
on the Earth. From Figure 1.1 we see that these include; (1) the 
transmitting antenna, hereafter called the spacetenna . (2) the receiving 
antenna, hereafter called the rectenna , and (3) the pilot signal trans - 
mitter . The rectenna and pilot signal transmit ter are located on the 
Earth. The purpose of the spacetenna is to direct the high-power beam 
so that it comes into focus at the rectenna. Since the rectenna is tc 
be approximately 10 km in "diameter," the beamwidth of the high power 
beam must be extremely narrow, in fact, on the order of one-half 
minute of arc. The pilot signal, transmitted frcxn the center of 
the rectenna to the spacetenna, provides the signal needed at the 
SPS to focus and steer the power beam. As a by-product of this 
requirement a command link to the SPS is provided. 

As seen from Fig. 1.1 the SPS phase control system is faced 
with several key problems. They include: (1) path delay variations 

due to imperfect SPS circular orbits, (2) ionospheric effects, 

(3) initial beam forming, see Figs. 1.2 and 1.3, (4) beam pointinn, 
see Fig. 1.4, (5) beam safinq. Figs. 1.2 to 1.4, (6) high power 

amplifier phase noise effects, (7) interference (unintentional and 
intentional), etc. 
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Figure 1.1. Space Based Solar Power Satellite and Earth Based 
Energy Collection System Concept. 
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1.1 SPS Transmitting System Concept 

From the Sjitem engineering viewpoint, the SPS transmitting 
system which Incorporates retrodl recti vlty Is depicted In Figure 1,5. 
A central feature of the SPS Transmitting System Is 101,552 element 
retrodl recti ve active phased-array (spacetenna) [2], [3] of 1-km 
diameter designed to focus and point the phase-coherent microwave 
beam to a ground-based rectifying antenna (rectenna) which Is 

approximately 10 km In diameter. As seen from Figure 1.5 the SPS 
Transmission System consists of three major systems; (1) The 

Reference Phase Distribution System , (2) Th e Beamforming and 
Microwave Power Generating System ,and (3) The Solar Power to 
Electrical Power Conversion System . In what follows we shall be 
concerned with establishing the far-fleld performance associated 
with the Reference Phase Distribution System and the Beamforminq 
and Microwave Power Generating System. The major technique used 
In the performance evaluation Is that of analytical simulation , i.e., 
a combination of analysis and computer simulation. 

1.2 The Phase Control System Concept 

To achieve retrodl recti vlty, the microwave power transmission 
system requires an on-board Phase Control System (Figure 1.5) phase- 
locked or synchronized to the pilot signal transmitted from the 
center of the ground based rectenna. One major purpose of the 
Phase Control System is to coherently reconstruct the Instantaneous 
frequency and phase of the received pilot signal and use this to 
generate and distribute a set of constant phase conjugation refer ence 
signals . This set of signals Is used to advance the phase of the 
pilot signal received at each spacetenna element by an amount equal 
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to the phase delay accumulated on the uplink. By a constant phase 
conjugation reference signal is meant a signal whose instantaneous 
phase characteristic, at any spatial point (x^.y^.z^) on the 
spacetenna, is perfectly correlated with the phase characteristic 
distributed from the spacetenna center to any other point 
on the spacetenna. Such instantaneous "phase-matching" in time and 
space implies synchronization [4], i.e., bringing about and maintain- 
ing a synchronous state at all points on the spacetenna where 
phase conjugation is to take place. In essence, the generation of a 

set of constant phase conjugation reference signals from the received 
pilot signal is equivalent to transfering precise time and time 
interval to a set of remote points (subarray centers) located on 
the spacetenna. At these subarray centers we locate Phase Control 
Centers (PCCs). In Reference 1 we illustrated and discussed their 
location on the SPS. 

1 .3 Beamforming and Microwave Power Generating System Concept 
To form a coherent, high power microwave beam a cooperating 
Microwave Power Generating System is required, see Figure 1.5. One 
of its main func‘:ions is that of delivering large amounts of micro- 
wave power to the radiating subarrays. Refs. 4, 5. By taking advantage 
of the retrod i recti ve feature of the active phased-array the power 
generated by the high power amplifiers is automatically returned 
(transponded) to the direction from which the pilot signal came. 

Thus the SPS system directs the power generated by the Space-Based 
Energy Station to the Earth based rectenna. The retrodirective phase 
conjugate concept, currently selected as the reference system. 
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win be further elaborated upon in what follows. Needless to say, 
there are other approaches to beam forming and power transmission; 
viz., ground control techniques which are currently being investigated 
at LInCom. 

1.4 The Reference Phase Distribution System and the Microwave 


stem Interface 




As previously mentioned the phase conjugate reference signals 
generated by the Reference Phase Distribution System are coherently 

distributed to Phase Control Centers (PCCs) located on the SPS. 

From these PCCs, the signals are then distributed and used to drive 

theset of phase conjugate circuits associated with the Beam Forming 
and Power Generating System, Figure 1.5. These processors, which 
contain the phase controlled high-power amplifiers operating in unison, 
automatically form the high power beam by causing advancement in the 
phase of the power amplifier output signals. This advancement is 
such that the accumulated phase on the downlink signal equals the 
phase accumulated on the uplink signal. Accomplishing this " zero 
phase accumulation orocess" on the up and downlink signals requires 


that the uplink frequency of the pilot signal "appear" to equal the 
downlink carrier frequency; this creates the SPS Pilot Signal Design 
problem. 

2.0 THE SPS PILOT SIGNAL MODULATION FORMAT FOR A RETROFIRE ARRAY 


The major factors of concern in the design of the SPS pilot (see Fig 2.0) 
waveform include: (1) ionospheric effects, (2) beam squint, (3) isola- 

tion of the up and downlink waveforms, (4) phase noise associated with 
the high power amplifiers, (5) diplexer characteristics, (6) security, 

(7) multiple acess capability and anti-jamming margin, (8) privacy 
of a command link to the SPS, (9) radio frequency interference 
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Figure 2.0, SPS Pilot Signal Design Considerations 
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Including that from other signals which may be generated in the SPS 
network (multiple access), (10) bandwidth requirements, (11) oscillator 
Instability, and (12) phase ambiguities generated in the Phase Control 
System. Appendix I provides analytical support to the selection of 
SPS pilot signal. 

Figure 2.1 summarizes the single-frequency pilot tone/retrofire 
concept investigated in Ref. 8. The major disadvantages of this approach 
are also summarized in Fig. 2.1. To overcome some of the disadvantages 
a two-frequency pilot-tone/retrofire concept has been considered. 

This approach is summarized in Fig. 2.2 along with its major disadvantages. 
Appendix I provides an analysis of the ambiguity resolution problem 
associated with the two-frequency pilot tone concept. 

To overcome the uplink pilot signal interactions with the high 
energy downlink power bram a pilot waveform has been designed (Ref. 1) 
which isolates the uplink from the downlink while simultaneously 
providing for anti-jam protection (security and privacy) and multiple 
access operation of any reasonable number of SPSs. In addition, the 
pilot waveform design has the added attributes that it tends to 
minimize the deleterious effects (beam squint and jitter) produced 
by the ionosphere's index of refraction. The modulation characteris- 
tics of the reference system pilot waveform is Binary Phase-Shift Keying 
(BPSK) of the pilot carrier for comnand transmission to the SPS and 
split-phase (Manchester coded), direct sequence (DS) pseudo-noise 
modulation for uplink/downlink Isolation and anti-jam protection [4,7], 
Code-division multiple access (CDMA) modulation is used to provide the 
needed multiple access capability required in the operation of a 
network of SPS's, see Figure 2.3. The SPS Pilot waveform design is 
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summarized as NRZ/BPSK/B I -<> -OS/CDMA, Ref. 1. i.e., for simplicity the 
signal is a Bi-^i coded tone, see Fig. 2.3. 

The math model for the SPS pilot signal s(t,c) is characterized 
as the product of three complex random modulations, i.e., 

Sy(t,;) • Re{v^ d(t,;)c(t,;)exp[ju)Qt+0(c)3} (2.1) 

P ^ ^ 

Command Code Carrier Modulation 

Modulation Modulation 

Here P^ represents the average transmitted power. In the notatiioni'M. i; i,H’.\iJTV 

used here, (e.g. d(t,;)), ; simply indicates that the quantity 

involved should be viewed as being randcxn in some way. A functional 

diagram indicating the mechanization of the pilot signal transmitter 

is shown in Fig. 2.4. As illustrated the data clock and code clock 

are coherent so that the uplink operates in a data privacy format. 

If the switch at the output of the RF oscillator is moved to the 
alternate position shown, then the data clock, code clock and 
carrier are all coherent. This configuration would simplify the 

I 

implementation of pilot signal receivers located at the SPS; the 
concern could be transmitter, media and oscillator phase instabil- 
ities. The purpose of the spread spectrum (SS) code generator is 
several fold. First it provides link security, second it provides 
a multiple access capability for the operation of a network of SPSs, 

See Fig. 2.5, and third, the anti-jamming protection can be provided j 

if needed. Proper choice of this code modulation will also provide 
the needed isolation between the uplink and the downlink. The 

j 

expression for the SPS transmitted signal s(t,c), indicated in (1.1), i 
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would be the same even If the uplink signal does not employ the 
spra- 1,1 spectrum modulation c(t,;), I.e., c(t,c) ■ 1 for all time. 

We snail refer to the signal c{t,;) as the SPS spread spectr um (SS) 
code signal. To simplify the demodulation process required In the 
spacetenna receiver, the SPS spread spectrum codes (multiple access) 
will product phase modulation of the complex carrier exp[J(u»Qt+e(c))3, 
I.e., 

lc(t,c)ls 1 (2.2) 

A mathematical characterization of the code Is given by 


c(t.c) » 2 *n**c^^“"^c^ 

n 

where {a } ■ 1 for all n and represents the chip time. Such an 
expression models Direct Sequence Pseudo-Noise Code modulation [ 1] 
and contains no random parameters. The notation p^(t) Is used there 
to denote a square T. second pulse of unit amplitude centered at the 
time origin and hence It Is easily seen that the direct sequence (DS) 

signals sstisfles the constant power and constant envelope condition. 
The complex sequence must be agreed upon in advance by the 

pilot signal transmitter and the SPS pilot signal receiver. In fact. 
It has a status similar to that of a loused in a cryptographic 
system. Thus with knowledge o^ the appropriate sequence, demodulation 
is possible and vfithout knowledge of that sequence, demodulation is 
extremely nifficilt !t is of interest to make the DS code purely 
random with no mathematical structure; however, all systems have a 
finite memory constraint and thus all practical OS codes have some 
periodic structure. I.e., a^ “ a^^j^ for all n where L denotes the 
period of the SPS code sequence. 
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2.1 Power Spcctriwi of the SPS Pilot Signal 

The term spread spectrtmi comes from the fact that the spectral 
density of the key signal c(t,;) Is very imich wider than the spectral 
density of the data modulation d(t.c). This means that the spectral 
density of the pilot signal s(t,;} Is much wider than It would be In 
the absence of the key c(t,;). 

If we now let z(t,;) » d(t,c)c(t,;) then 

{*/?PJJ)’^s(t,;) ■ Re{ 2 (t,c)expCJ(wQt+e(;))]} (2.4) 

and under reasonable assumptions about the process d(t,c) and c(t,i;) 

It Is easy to show that the power spectral density of the amplitude 
process z(t,;) Is given by 

S^(f) - S^(f) * S^(f) (2.5) 

where S .(f) and S.(f) are the power spectral densities of the data 

y C 

modulation and OS code modulation respectively. Here denotes 

the convolution operation. We now discuss the power spectral densities 
of the data modulation and DS code modulation respectively. Here 

'•*" denotes the operJtion convolution . We now discuss the power 

spectral densities for Sjj(f) and S^(f). 

2.2 Power Spectral Density of the Data Source 

A convenient representation for the SPS uplink command modulation 

is 

d(t,c) ■ ^ 

n«-« 

where c Is a random pulse epoch uniformly distributed over the ^ 
duration T^^ seconds. Additionally, assumed to be a 






om 






j 


sequence of In^pendent, equelly probable conmand delta bits taking 
on values plus one or minus one. The uplink command data rate Is 
therefore 


R 


b 


J- bits/sec 
'b 


(2.7) 


Using the Wiener-Khintchin theormn, It Is easy to show that [7] 




(sin irfTjj)^ 


( 2 . 8 ) 


This power spectral density Is Indicated In Figure 2.6. Notice that 
the main lobe occupies a bandwidth of approximately B • 2 / 1 ^ Hz. 

2.3 Power Spectral Density of the OS Code Modulation 

If the SPS code modulation c(t,;) Is to be used it should possess 
the following features: 

(1) Wide power spectral density in order to achieve anti-jam 
(AJ) margins if required. 

(2) Ito line spectrum cofl^onents to prevent jamming by 
narrowband Interference. 

(3) Spike-1 Ike autocorrelation function for accurate code 
synchronization. 

(4) Non-predictability for use against intelligent jammers 
and to maintain cwnplete security and privacy on the 
uplink modulation. 

(b) Low cross-correlations with other possible DS codes to 
render the multiple access capability reouired in the 
SPS network. 

(6) Possess spectral characteristics such that the SPS uplink 
and downlink signals do not lit. in the same frequency. 



1 

i 

i 

j 
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(7) Possess spectral characteristics so as to isolate the 
SPS uplink signal from the downlink SPS high power 
amplifier output phase noise characteristics. 

(8) Possess spectral characteristics such that the effects 
produced by the ionosphere delay on the upper sideband 
frequency components averaged with those of the lower 
sideband frequency components equals the average ionospheric 
effects produced on the carrier component at frequency fg Hz. 
This minimizes the effect of the varying index of refraction 
produced in the ionosphere on the downlink power beam. 


Particular types of random key sequences most frequency used 
for constructing DS code keys are the so-called pseudo-noise (PN) 

sequences possessing a period of 2*"-l code chips. The discrete 
periodic correlation function, say R(.(t), of a PN sequence {a^} is 

defined as 

2*--l a. ©a. 


R,(t) • Y. " 
n=l 


n-T 


where © denotes mod 2 arithmetic in the exponent. Using the properties 
of PN sequences [7,8] the correlation function becomes 


R,(t) 


2*--l 


-1 


T = 0 
T 0 


Since the sum is over one period of a periodic sequence, the result 
is independent of the starting state of the PN sequence. The code 


,N 


period is actually (2 -1 )T^ seconds in duration. 
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In order to satisfy requirements (6), (7), and (8) in the above 
list an additional modulation must be applied to the generated code, 
viz., that of bi -phase (Bi-cj)) modulation of the code chip pulses. 

Ref. 14. This is sometimes referred to as Manchester coding and 
is most frequently applied to the data source pulses for the purposes 
of guaranteeing a certain data transition density in the data modula- 
tion d(t,c) for purposes of clock recovery. Applying this modulation 
to each code chip pulse results in the modified code generated output 
waveform of 

c(t,c) = ^ a^(c)p^(t-nT^)Pj^(t-nT^) (2.9) 


where P^(t) constitutes the Bi-ij) pulse waveform 


Pjt) = u(t) - 2u(t-T^/2) + u(t-T^) (2.10) 

and u(t) is the unit step waveform. Using the Wiener-Khintchin 
theorem it is easy to show, under suitable assumptions, that the 
power spectral density of the Bi -phase code key is [7] 


5,(f) = T 


c 


sin^(iTfT^/2) 

(^fT^/2)^ 


( 2 . 11 ) 


where T is the code chip duration. This power density is illustrated 
c 

in Figure 2.2 for the case where the code processing gain 

PG = 1^/T^ (2.12) 

is much greater than one. Now approximately 90°$ of the area under 

2 

the (sin x/x) curve lies in the central peak between the f'‘rst 


zeros of the data modulation spectrum and the first side lobes on 
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both sides contain another S% of the area. If PG >> 1, then the 
spectrum of d(t,c) is much narrower than that of c(t,c); hence, 
S^(f) behaves approximately like a delta function when convolved 
with S (f). Thus from (2.8), (2.11), and (2.5) we have 


s^(f) 


sin^(irfT^/2) 


(^fT^,/2) 


T 


(2.13) 


while the power spectral density of the transmitted SPS pilot signal. 


as found from (2.4) and (2.13), is 


Sg(f) = 


^u^c 


sin^7r(f-fQ)T^/23 sin^[^(f+fQ)T^./2] 


(2.14) 


This spectral shape is illustrated in Figure 2.7. Notice from (2.14) 

that as the chip time T^ becomes smaller the spectral null around f^ 

broadens and the overall power spectral broadens (features (1) and 

(2) in Section 2.2) and decreases in amplitude. Thus the code chip 
duration T^ serves as the parameter which can be used to control 

the degree of uplink/downlink isolation, see feature (6) in 

Section 2.2. Variation of T^ can also be used to tradeoff uplink 

bandwidth requirements with diplexer and other receiver filter 

requirements. 

If the switch appearing at the output of the code generator 
were switched to the alternate position then the uplink is no 
longer secure nor would it possess the multiple access capability; 
however, the pilot signal still retains the required spectral and 
sideband averaging properties needed to isolate the up and downlink 
signal. The transmitter output power spectrum for this configura- 
tion is sketched in Fig. 2.8 for the case where the transmitted 
carrier is not modulated with the spread spectrum code. The split- 
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phase format of the PN (pseudo-noise sequence results in a spectral 
distribution of energy which peaks either side of the carrier at 
approximately 3/ 47^, where R Is the PN code chip rate. Proper 
selection of will thus allow RF diplexing at the power transponder 
receiver to provide isolation between the desired pilot signal and 
the unwanted power signal which Is present at the pilot carrier 
frequency of 2,450 MHz. 

2.4 Code Division Multiple Access (CDflA) Modulation 

The spread spectrum code c(t,;) possess an Inherent multiple 
access capability which will be needed in the operation of a network 
of SPSs. For example, each SPS can be assigned a pseudo-random 
(PN) code (or one from the so-called Gold family) each of which 
Is orthogonal (or near orthogonal) with respect to each of tie 
other assigned SPS codes. In fact, a single PN code generator* 
design can be used for all SPSs with differing logic combinat ons 
used to derive the desired set of orthogonal codes. 

CDMA techniques therefore offer an attractive multiple access 
approach and some anti-jam protection if desired. They require 
little or no additional equipment complexity in their generation 
and the bandwidth occupancy requirement Is equal to that required 
for a single SPS. The geometric dispersions of the SPSs should not 
be overly severe to Impose restrictive constraints on CDMA 
performance. 

3.0 REFERENCE SYSTEM S PS PILOT WAVEFORM SUMMARY 

The reference system SPS pilot waveform utilizes: (1) NRZ command 

modulation, (2) split phase, direct sequence pseudo-noise or spread 


spectrum modulation, BI- 41 -DS. This combined data-code modulation is 
used to bi-phase modulate (BPSK) the RF carrier. Multiple access in 
the SPS network is to be achieved via code division multiple access 
techniques (COMA). Thus the baseline SPS pilot waveform is characterized 
via four modulation products summarized by the symbols: 


Command 
Modulation 


NRZ / BPSK / BI-<>-DS / CDMA 

t: 


T 

RF Carrier 
Modulation 


Spread ~ ^Multiple Access 
Spectrum Modulation iModulation 


The coimand data rate, R^^ * l/T^^ bits/sec, nor the code chip 
rate, R = 1 /T chips per second, is not specified in this report; 
however, LinCom has j'’:;t developed preliminary chipping rates. 

The baseline SPS up and downlink signals are designed to occupy a 
band of frequencies centered at fg = 2450 MHz; thus the bandwidth 
requirement is currently under study at LinCom. 

4.0 R EFERENCE PHASE CONTROL SYSTEM 

As previously noted, the phase control system (see Fiq. 1.2) 
must accomplish three major functions. These include: (1) power 

beam forming, (2) power beam pointing, and (3) power beam safing. 

The first involves the generation of a highly directional narrow 
beam microwave signal. This is accomplished at the spacetenna by 
coherently phasing each radiating power module (subarray) in the 
one kilometer diameter spacetenna to produce the desired broadside 
spacetenna far-field pattern. In addition, in order to improve 
power transfer efficiency, the power transmitted over the spacetenna 
"spelling" aperture is tapered by approximately 10 dB from the 
spacetenna center to its outer edge. This power taper effects the 


beam shape by concentrating more of the microwave energy Into the 
center portion of the main beam by widening the main lobe of the 
antenna pattern. Later on we present some average power patterns 
resulting from the computer simulation of the phase control system. 

Once a coherent beam Is achieved, the center of the beam must 

be precisely pointed at the Earth-based rectenna. Since the 

mechanical pointing accuracy of the spacetenna and the Individual 

subarrays within the spacetenna Is expected to be on the order of 

one to three minutes of arc, the potential miss distance on the 

Earth's surface from the geosynchronous orbit would be on the 

order of 10 to 30 kilometers which would completely miss the 

10-k11ometer diameter rectenna site. To compensate for this 

pointing Inaccuracy, the phase control system of the SPS must be 

capable of adjusting the phase of each radiating element power 

module (subarray) to shift the power beam center without degrading 
the beam shape. 

The final function to be provided by the phase control system 
Is power beam safing. This aspect is characteristic of the phase 
control process selected since the power beam intensity becomes 
diffused If phasing of the high power amplifiers Is uniformly 
random. As we i’.hall see complete loss of phasing results In increased 
sidelobe levels, the resultant Increase remains well below the USA 
and USSR standards. The other aspect of dephasing has to do with 
power robbing by Intentional or unintentional jamming. By utilizing 
the concept of a retrodirective phased-array incorporating the pilot 
signal design discussed in the previous section, and perhaps ground 
sensors for additional safing considerations, the phase control 
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system sunmarized In the next section meets the general functional 
requirements. 

4.1 Phase Control Via Retrodi recti vlty 

The phase control system concept to be simulated was presented 

In detail In Ref. 1; Its major features will be summarized In this 

section. Based upon earlier study efforts (Refs. 1,8), a phase 

control system concept has been proposed which partitions the system 

Into three major levels. Figure 4.1 demonstrates the partitioning 

and represents an expanded version of Fig. ’.2. The first level In 

Fig. 4.1 consists of a reference phase distribution system Implemented 

In the form of phase distribution tree structure. The major purpose 

of the tree structure Is to electronically compensate for the phase 

shift due to the transmission path lengths from the center of the 

spacetenna to each phase control center (PCC) located In each 

subarray. The second level Is the Beam Steering and Microwave 
Power Generation System which houses the SPS Power Transponders. 

This transponder consists of a set of phase conjugation multipliers 

driven by the reference phase distribution system output and the 

output of a pilot spread spectrum receiver (SS RCVR) which 

accepts the received pilot via a diplexer connected to a separate 

receive horn or the subarray Itself. The output of the phase 

conjugation circuits serve as Inputs to the third level of the 

phase control system. The third level of phase control Is associated 

with maintaining an equal and constant phase shift through the micro- 

wave power amplifier devices while minimizing the associated phase 

noise effects (SPS RFI potential; on the generated power beam. This 

is accomplished by providing a phase-locked loop around each high 

power amplifier. The amplifier phase Is therefore stabilized by 
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phase locking to the conjugated phase reference provided at the 
output of the phase conjugation multipliers. 

4.2 Optimal Topological Design of a Reference Phase Distribution 


Feed Structure 

Based upon an earlier UnCom study (Ref. 8 ) various techniques 

to the problem of generating, maintaining and distrlbutina a coherent 

reference signal over a large area were considered. These techniques 

were modeled and mathematically analyzed with respect to their 

ability to minimize: phase-error build up, beam diffusion, beam 

pointing Jitter, cable length and maximize power transfer efficiency. 

To accomplish this distribution a tree structure has been selected; 

In fact. Appendix II, presents preliminary results related to the 

optimal design of a survivable feed structure for the reference 

phase distribution system. As seen from Fig. 4.1, the reference 
distribution tree consists of three major parameters and one functional 

building block. The parameters of the tree to be optimized include: 

(1) the number of levels in the tree, (2) the number of branches per 

level and (3) the interconnecting cable length. The functional 

building block consists of the Master Slave Returnable Timing 

System (MSRTS) units which is briefly discussed in what follows. 

4.2.1 The Master Slave Returnable Timing System (MSRTS) 
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a portion of the reference signal at the received location Is amplified 
and coupled back Into the phase distribution cable via a circulator. 
This return signal component Is then used for phase locking to the 
original reference Input of the Phase Control Center (Prc). Thus 
the phase of the VCO in the PCC Is advanced by an appropriate amount 
so that the phae of slave signal e 2 (t) equals the phase of the 
reference signal e^(t) transmitted via the "cabl^," see Appendix II. 
Since the same frequency is transmitted both directions over the 
cable, phase error build up due to frequency dispersion on the cable 
Is elimlnateu. When phase lock Is achieved in the PCC, the resultant 
phase at point B In Fig. 4.2 equals the phase at point A even if the 
effective cable delay varies. If the connection cable Is cut to 

within plus or minus one-fourth of an Integer number of wavelengths 
at the distribution signal frequency, the phase can be distributed 

accurately over various lengths of cable without phase ambiguity. 

For example, for a frequency of 490 i'friz, the cable lengths must be 
cut to within plus or minus fifteen centimeters to avoid ambiouities. 
4.2.2 Phase Transfer Between the Phase Control Centers 


The MSRTS discussed In the previous section is u«ed as bu^ldinq 
blocks to transfer the phase of the SS pilot reference receiver 
output to various points on the spaceteoiia where phase conjuoation 
can take place. Since each building block in the MSRTS contains 
a PCC, see Figures 4,1 and 4.2, onemust specify the pCC locations 
and ultimately how the phase conjugator reference signals are to be 
supplied as inputs to the Power Amplifier Conjugation Electronics 
which feeds the radiating elements. There are three basic confioura- 


*By "cable" is meant any suitable propaoation medium connect i.ig the two 
desired points. 
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tions which use the MSRTS concept to transfer the phase between PCCs. 
This includes: (1) Series Fed PCCs, (2) Parallel of Tree Fed Pees, 

(3) Hybrids of the first two, and (4) Attachment of PA’s at each 
tree level. In what follows we discuss the parallel or tree fed 
structure which has been chosen for use in the baseline. 

4. 2. 2.1 Tree Structure Used to Transfer Phase Between PCCs 
Figure 4.3 illustrates a typical tree structure connecting 

the PCCs wr:' .1 transfer and maintain the reference signal phase 

character! St’ c over the spacetenna aperture. As can be seen from 

Figure 4.3, the pilot reference signal (R) is transferred to four 

PCCs appearing at level 1 in the tree. Each PCC in level 1 is 

connected to n PCCs located at the next level in the hierarchy. 

Each PCC at this level feeds n PCCs located at the next level in 
the hierarchy. This process is repeated until the required number 

of PCCs is fed. 

In general, a symmetrical tree is defined by two parameters, 
viz., the number of levels, say L, and the number, n, of branches 
per level. For this case there would be generated n*" tree nodes 
(PCCs). In such a case, the tree structure could be used to provide 
n*" phase conjugator reference signals. Here the parameter "n" 
defines the number of output terminals associated with power dividers 
required in the implementation. 

4 . 2 . 2 . 2 Tree Structure Generated Using the MSRTS Configuration 
Figure 4.4 illustrates how the phase characteristic (derived 

in the reference SS receiver) can be transferred using the MSRTS. 

Here we have selected the MSRTS which uses a single frequency and 
circulator shown in Figure 4.2. The 490 MHz signal provided by the 
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Costas loop in the reference receiver of Figure 4,1 is frequency 
multiplied by two and used to drive the PCCs located at the first 
level in the tree, see Figure 4 .4. As seen each MSRTS feeds a 
power splitter and a times two multiplier at the next tree level. 

Tnis partitioning continues until the tree structure is completed. 

At the final stages of the tree the set of phase conjugator signals 
required for driving the conjugation circuits in Fig. 4 .1 can be • 
generated. Since there are 101,552 power amplifiers and 101,552 
radiating/receiving elements are specified then the parameters n 
and L must be selected so as to best "match" these characteristics 
to the tree. Once this "match" has been determined one must spec ify 
the layout on the SPS. Performance tradeoffs between the parameters 
are presented later on in this report. 

4.2.3 The Reference Phase Distribution System Functional Diagram 
Figure 4.3 demonstrates the first three levels in the layout 
of the reference phase distribution system. Here we illustrate 
reception of the pilot signal at the center of the spacetenna. This 
signal passes through the diplexer (D) to the reference spread 
spectrum receiver (R). The SS receiver reproduces the carrier 
component in the received pilot signal. The frequency of this signal 
is doubled to produce the 980 MHz reference. The reference is sub- 
sequently power divided and supplied to four PCCs via the MSRTS, 
see Fig. 4.5. The MSRTS process is replicated until the appropriate 
number of levels in the tree are developed. The outputs of the tree 
then serve as phase conjugator reference signal needed for beam 
steering and for stabilization of the high power amplifier phase 
characteristics by means of a phase-locked loop. 



Figure 4.4. Phase Control System Configuration Near Spacetenna Center 


































rhree Levels of the Referen'-e Phase Distribution System 
Illustrated in all Componer s. 
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The number of levels In the tree, L, as well as the number ot 
branches per level, k, represents parameters to be optimized from 
a far*f1eld performance point of view. As we shall see, the rms 
pointing error and the main beam gain reduction are affected by 
these design and layout of the tree structure. Thus the far-field 
performance for various tree structure designs. The redundancy 
requirements and possible other mechanical layouts are under investi- 
gation to insure .that the reference phase distribution network will 
exhibit high survivability and reliability characteristics, see 
Appendix I. 

3.0 REFERENCE SYSTEM SPS POWER TRANSPONDER 

In addition to distributing the constant phase reference signal 

over the spacetenna as discussed in the last section, a method for 
recovering the phase of the received pilot signal is required. 

Figure 5.1 represents the functional diagram of the SPS power 
transponder. This includes the pilot signal receiver, phase 
conjugation electronics and the high power amplifier phase control 
subsystem. 

In the mechanization of the SPS power transponders, two 
receiver "types" will be required; however, most of the hardware 
will be common between two receivers. One receiver, the Pilot Spread 
Spectrum Receiver , is located at the center of the spacetenna or the 
reference subarray. It serves two major functions: (1) Acquires 

the SS code, the carrier and demodulates the command signal, (2) 
provides the main input signal to the Reference Phase Distribution 
System, see Figure 5.1. 
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The second receiver "type" will be located in the Beam Forming 
and Microwave Power Generating System, see Figures.^. Its main 
purpose is to phase conjugate the received pilot signal and trans- 
pond power via the spacetenna element, j * 1 ,2,. . . ,101 ,552. We 
I now discuss the functional diagram indicating the mechanization of 

the SS power transponder and discuss its operation. 

^ From Fig. 5.1 we note that each SS receiver element must be 

capable of despreading and demodulating the received pilot signal 
based upon the gain provided by a single element of the SPS 
spacetenna. The receiver consists of several major subsystems. 

These include, (see Fig 5.1 and 5.2): (1) the RF Front End centered 
j at 2450 MHz, (2) The SS Code (PN) Sync Subsystem (PNSS), (3) The 

AGC Subsystems, (4) The Carrier Sync Subsystem (CSS), (5) The 

Carrier Lock Detection Subsystem, (6) Carrier Sync Acquisition 
Subsystem, (7) Symbol Synchronization Subsystem, (8) The Viterbi 

Decoder (if the uplink employs convolutional coding), and (9) The 

Ambiguity Resolving Subsystem. 

The philosophy used in partitioning the receiver is summarized 
below: 

•Rapid acquisition in the Carrier Sync System 

I •Despreading Process Independent of Carrier Acquisition, 

Tracking and the Demodulation Process 

The carrier, with nominal frequency of 2450 MHz, is first processed 
via an RF filter. The bandwidth of this filter must be sufficiently 
wide to pass the PN chips and the roll-off must be sufficient to meet 
the desired rejection requirements. 
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Power Transponder 

















The signal level Into the first IF mixer Is held constant by the 
noncoherent AGC . The first LO Is selected to run at 1960 MHz; 
therefore, the first IF frequency, at zero Doppler, Is 490 flHz. The 
output of the first IF mixer Is further filtered by the first IF filter 
whose bandwidth Is sufficiently wide to pass the PN chips. The cascaded 
frequency response of the RF filter and the first IF filter are collect- 
ively designed to meet the desired front end rejection requirements. 

The IF filter output serves as the Input to the PN synchronization 
system. The PN synchronization system (PNSS) of Figure 5.1 incor- 
porates a noncoherent -t-dlther PN acquisition and tracking desian. 

The details associated with the particular algorithms selected for 
the code acquisition will not be discussed here. 

The PNSS is configured to give rapid acquisition of the chosen 

PN code in noise and in the presence of worst-case Doppler frequency 
shifts. The PNSS also provides highly reliable code loop tracking 

performance under a chosen minimum Input signal-to-noise level without 
severely degrading the symbol synchronizer, carrier tracking loop and 
data detection processes. 

As shown in Figure 5.1, the output of the second IF mixer 
(despreader) is processed via two paths. One path is selected for 
signal demodulation and carrier recovery while the other for Pf! code 
acquisition and tracking. Based upon hardware considerations a 
T-dither early/late gate code tracking loop has been chosen. This 
loop generates a code error signal proportional to the phase difference 
between the locally generated PN code and the received signal code. The 
code error signal generated at IF is filtered in the bandpass filter/ 
envelope detector arrangement shown, A loop filter integrates the 


code phase error signal and drives a VCO to null the error signal. This 
assures precise code alignment at the on-time correlator or PN despread- 
ing mixer. The PN despreading mixer output contains the data signal 
without the code and Is then processed as a normal BPSK suppressed 
carrier signal by means of the CSS, see Figure 5.1. 

The PNSS noncoherent code acquisition design point slgnal-to- 
nolse ratio has not been determined; however, a code phase synchronization 
search would be performed over the entire doppler-expanded code length 
by stepping the local PN code generator at half-chip Intervals, Integrat- 
ing the envelope detector's output for a dwell time of T^ seconds, and 
then making a sync decision based upon a threshold comparison of this 
Integrated output. The threshold Is established by means of a noise 
estimate produced via an appropriately chosen algorithm. If the decision 
is not In favor of the particular code phase being searched, the code- 
tracking loop filter Is quenched and the search proceeds to the next 
code phase position. The PN code acquisition stage of synchronization 
ends with the two codes (transmitted and local) most probably 
aligned to within an average of one-quarter of a chip. The second- 
order early-late, i-dither, code tracking loop will then pull the 
system into final alignment. 

The code error signal is recovered at the output of the bandpass 
filter/envelope detector shown in Figure 5.1. The error signal 
polarity Is toggled in synchronism with the early-late switchinn of 
the local PN code thereby creating a voltage proportional to the 
phase difference between the locally generated and received PN wave- 
forms. This voltage is filtered by the loop filter and used to 
drive the VCO so as to null the error signal. Once the code loop 
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locks, the loop filter positions the code loop VCO frequency so 
that the phase error between the local and received codes approach 
zero. 

The arriving NRZ/BPSK/B I -OS/COMA signal Is despread prior to 
filtering by the second IF filter. This gives rise to an ordinary 
BPSK signal when the r-dlther loop Is locked. The output from the 


second IF filter goes to the output of the carrier recovery circuit. 
The bandwidth of the second IF filter Is chosen to be wide with 
respect to the data rate. This avoids creating Intersymbol 
Interference In the data stream thereby minimizing command data 
detection errors, avoiding false lock problems In the suppressed 
carrier tracking loop as well as minimizing phase jitter on the 
reference signal which drives the Phase Control System. 

A Costas (I/Q) loop configuration 1s chosen for carrier 
acquisition, tracking and data demodulation. This configuration 
<as chosen because It was determined to be optimum when all 
considerations. Including the ability to square perfectly over 
temeprature and signal level, are traded against lock detection 
and synchronization monitoring. We note that the bandwidth of the 
Costas loop serves to control the time dependent phase jitter fed 
into the Phase Control System. In order to minimize this effect 
the bandwidth of the Costas must be as small as oscillator and 
media instabilities will allow. 

Noncoherent AGC is derived from the sum "I" and plus "g" 
Channels appearing In the arms of the Costas Loop. In addition, 
lock detection for the carrier circuit Is accomplished by usina 
the difference between the "I" and "0" channels or' the Costas amis. 


JUinCc 


om 


- 46 - 


The noncoherent AGCg controls the receiver gain (prior to phase 
detection) with the signal plus noise level appearing In the 
outputs of the Costas loop arm. This feature Is used to control 
the loop bandwidth and damping factor during acquisition and track- 
ing. The bandwidth of the arm filters are chosen to minimize the 
so-called squaring loss and their roll-off characteristics selected 
to avoid the false lock problem. 

The loop filter receives the signals from the phase detector 
(third multiplier) and supplies an error signal to the VCO which 
control the local frequency. The loop filter sets the trac^'lnq 
loop bandwidth and damping factor. An A6C voltage, proportional 
to the Incoming signal plus noise power, Is low pass filtered and 

amplified to drive variable gain elements In the IF amplifiers. 

(These amplifiers are Included here in the first and second IF filters 

for simplicity.) Additional Integrate and dump circuits and threshold 

detector circuits control the sweep and the lock Indicator needed 

for loop supervisory control. It Is to be noted that the carrier 

sweep Is not activated until PNSS Is synchronized. 

The carrier lock detector circuit monitors the Integrate and 
dunp voltage formed at the discrete points In time by differencing 
the squares of the Inphase and quadrature arm outputs. The output 
of the Integrate and dump circuit Is compared to a fixed threshold 
level. If the I and 0 output exceeds the fixed threshold level a 
preliminary Indication of lock Is given and the search sweep Is 
disabled. When lock Is verified the loop bandwidth could be reduced. 
If necessary, to further reduce loop jitter. In such a case the 
loop bandwidth would be optimized for both acquisition and tracking. 
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Data extraction can b<j derived in two different ways. The 

most convenient way Is to extract the data from the output of the 

Q channel and apply the symbol stream to the symbol sync system 

as Shown In Figure 5.1. The symbol synchronizer extracts the NPZ 

baseband waveform and applies the symbol (soft decision) voltages 

to the Vlterbi decoder (If coding Is used) for further processing. 

The 180 degree phase ambiguity In the Costas loop Is of no consequence 

If \ convolutional code (asswning the link Is coded fo* error control) 

Is selected that Is not sensitive to code polarity. On the other 

hand* If a code Is selected that Is sensitive to code polarity 

(transparent vs nontransparent code) then the decoder can be used 

to resolve the 180® ambiguity In the Costas loop. Whether coding 

is used or not, proper operation of the phase conjugation circuitry 
cannot be achieved unless this ambiguity is resolved. A method for 

achieving this Is provided In Appendix I. 

An alternate approach to extraction of the commjnd data Is 
to use a wideband phase detector and reflect the data to baseband 
by ucing the despread signal which appears at the input to the 
second IF filter. However, the differential phase shift accumulated 
between the second IF filter input and the input to the wideband phase 
detector created by the loop must be controlled if this abroach is 
used. It is not recommended here. 

The 490 fWz reference appearing at the pilot receiver VCO output, 
see Figure 5.1, serves as the only input to the Phase Control System. 
This phase characteristic is then distributed over the aperture o^ the 
spacetenna. The output of the Phase Control System consists of 101,552 
980 f'Wr constant phase reference signals which are used to conjuqtste 
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the incoming signal and to stabilize the high power amplifier (HPA) 

outputs. This is acheived by placement of an automatic phase control 

around each HPA, see Figure 5.1. This will be discussed further later. 

Overall transponder performance is affected by various system 

interactions. These include the AGCs, PNSS, CSS, SSS, and the power 

power amplifier phase control system. Their performance will not be 

further discussed here. A time domain simulation of the SPS receiver 

elements is currently available at LinCom for operation in both the 

acquisition and tracking mode. This simulation is currently being 

extended to include other elements of the Power Trandponder, viz., 

the Phase Control System and the Power Amplifier Phase Control System 

and associated phase conjugation circuitry. Figure 5.3 provides a cumnary of 

the "high level" electronic subsystems required in the implementation of the 
SPS Reference Phase Distribution System shown in Fig. 4.1. 


5.1 Conjugation in Groups of Power Modules 

The reference system described above needs a conjugator at each 
power amplifier and hence needs 101552 conjugators. An alternate 
way of doingthe same thing is to put a single conjugator at the 
center of constant size subarray (10 mxlO m), conjugate the received 
pilot phase at that point and then sending this conjugated phase to all 
the power amplifiers in the subarray via MSRTS. This arrangement 
simplifies the hardware implementation. If there are no mechanical 
errors, i.e., jitters due to random errors are negligible then this 
system works just as well as the baseline but this system is very sensitive 
to the random errors and tilts of one subarray. This is described in 
greater detail in Section 10 of this report. Figure 5.4 depicts the 
system concept and functional diagram which serves to perform the 
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(1) NUMBER OF SS RECEIVERS 

•COSTAS LOOPS 
•DESPREADERS 


101,553 

101,553 


(2) NUMBER OF DIPLEXERS 101,552 

(3) NUMBER OF POWER MODULES 101,552 

(4) NUMBER OF PHASE CONJUGATOR 

MULTIPLIERS 101,552 

(5) NUMBER OF 4-WAY POWER SPLITTERS 40,960 

(6) NUMBER OF MSRTSs 22,000 

(7) APPROXIMATE "CABLE LENGTH" 

REQUIRED* 120 miles 


*SEE FOOTNOTE ON PAGE 34. 


Figure 5.3. Summary of the High-Level Electronic Subsystems 

Required in the Implementation of Phase Distribution 
System, the Beamforming and Microwave Power 
Generation System. 
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conjugation at the center of each constant sized subarray. This 
simplification should be compared with the reference system concept 


in Fig. 4.1 . 

6.0 SPS PHASE CONTROL SYSTEM PERFORMANCE EVALUATION VIA THE METHOD 

or sTATisTicALTrrrmiA Theory 


Because of the complicated nature of the problem of evaluating 
performance of the SPS phase control system and because of the 
multiplicity and interaction of the problems as they relate to sub- 
system interfaces, the methods of analytic and computer simulation 
(analytical simulation) have been combined to yield performance of 
the SPS system. Figure 6.1 illustrates the approach which has been 
taken in order to develop the performance of the phase control 
system described in the previous sections of this report. As 
seen from this figure various inputs in the form of requirements , 
technology status, constraints, and other NASA/DOE contracts have 
been used in the problem formulation. From this formulation it has 
been convenient to break the analysis of the phase control system 
into three major problem areas. Outputs from these analytical 
efforts have been used to provide direction for development of a 
computer simulation (called SolarSim) whose major output is the 
performance of the SPS retrofire phase control system. The SolarSim 
outputs (to be discussed later) are then used to make System 
Engineering tradeoffs and perhaps support modifications in the 
baseline system. 

6 . 1 SPS Phase Control System Performance Measures 

In this section of the report we shall discuss the pertinent 
performance measures associated with evaluating the far-field 
characteristics and performance of the SPS antenna when its current 
excitation is driven by random fluctuations arising in the Phase 
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Control System, The randomness produces errors in the amp! itude 

and phase components of its element feed currents. 

The major sources of errors are due to: (1) ionosphere (propaga- 
tion medium), thermal , PA phase noise, etc,, (2) subarray defects 
and fabrication error, (3) phase error bui''dup in the phase control 
systen and (4) motion induced irregularities and thermal stresses. 

In conjugation with the possible error sources responsible for the 
random fluctuations, one must distinguish between the statistics 
taken over the ensemble of similar antennas and the statistics in 
time from an individual SPS antenna. 

In the first case, one investigates the scatter of the 
parameters of antennas of the same type with respect to the ensemble 
of antennas. Such a scatter may be caused by inaccuracies in the 
production of the feed structure, by the inhoniogeneity of the 
material of which the SPS antenna is made, etc. 

In the second case, one investigates the scatter of parameters 
of one and the same spacetenna in time. This scatter may be caused 
by instabilities of the feeder elements of the spacetenna. random 
changes of the spacetenna surface due to solar pressure or temper- 
ature changes, changes iii the parameters of the ionosphere through 
which the pilot signal has propagated, etc. 

The problem of analyzing the SPS spacetenna performance as 
effected by the Phase Control System is divided into two partial 
problems » viz., the internal and the external problem. The former 
problem deals with modeling of the statistics of the amplitude 
and phase excitation currents in accordance with the implementation 
of the phase control system. This problen is approximation methods, 
viz., modeling of the actual mechanisms responsible for error 


sources. On the other hand, the external problem consists of 
characterizing the statistics of the far-field using the solution 

to the internal problem. 

Using the statistical characterization of the far-field a 
number of important characteristics of the transmitted power beam 
can be studied. They include the spacetenna pattern and such 
associated parameters as gain, beamwidth, sidelobe level, the phase 
and polarization patterns. In this study we have concentrated on the 
evaluation of the mean (average) spacetenna far-field characteristics 
as a function of the phase amplitude error build-up in the phase 
control system. The approach has been to develop an analytical 
simulation capable of predicting the performance of the phase control 
system. By analytical simulation is meant the combination of usinn 
analysis and computer simulation to develop a computer nrogram, called 
SolarSim, to develop the far-field performance characteristics . 

6.1.1 Phase Control System Performance via S olarSim 

In order to assess the phase control system performance, it has 
been nv^cessary to put the multidimensional phase control problem 
under computer control via the development of the analytical simulation 
called SolarSim. SolarSim, see Fig. 5.2, is an SPS system engineer- 
ing computer simulation design'd to reduce the entropv of the SPS 
phase control system problem. A description of the SolarSim software 
is given in Appendix 1. Because of the inability to implement a full 
scale model of the SPS phase control system to demonstrate perfonn- 
ance, it appears that the only way to obtain performance of a full size 
system is to simulate a full-scale model. Thus the implementation 
of a reduced size model can be used to verify its performance 
concurrently with the program. Once the program is fully developed 
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Figure 6.2. SPS System Engineering Computer Simulati 
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In this section we shall define the various performance char- 
acterization of the SPS spacetenna in terms of the far-field pattern. 
The geometric center of the spacetenna coordinates of interest are 
illustrated in Fig. 7.1. The particular realization of the antenna 
pattern for a retrodi recti ve antenna can be written as 


jKj[ ( sin 6 cos 4i - sin 0 ^ cos 
k 


+ (sin 6 sin <li - sin 9 sin )y. ]{7.1) 

r r kn 


With 


sin 0 cos 4 = u sin e cos 4 = u 

r r r 


sin 0 sin “ V 


sin cos 4 ' V 
r '^r r 


This can be rewritten as 
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where represents the total phase error build up at the (k ,0 


(7.2) 

th 


radiating element having coordinates (Xj^.y^^) and is the magnitude 

,th 


of current associated with (k,ii)— radiating element. This amplitude 
factor may indeed be random. Notice that in this formula the effect 


of phasing the tree structure is not explicit hence we have 





f(u,v) - I 
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( 7 . 3 ) 


where t|»_ _ signifies the random phase generated by the distribution 

in^ • • • 

tree by going through the branch at the 1— level 1 ■ 1,...,L. 

The power pattern Is then defined as f(u,v)f*(u ,v) » |f(u,v) ,1 .e. , 
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while the average power pattern becomes 
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Here we have assun.. ' that the current amplitudes and phases are 
Independent random processes. The effect of the randomness Introduced 
by the phasing tree is clearly seen In this equation. Note that the 
difference of phases between elements Is more important than the 
absolute phase of the element. Also note that the currents are 
only the amplitudes. From this average power pattern one can observe 
the beamwidth, sidelobe levels, shifting of the nulls, variation of 
sidelobe levels. 



- 60 - 



n 







om 


'oLin 




where 



m I . 




.m 


L 





• (<l> / 

' . . .iTi|^ m-j . . 

\*^v’^v ' ’®u C^respectively , but with x and y 

interchanged. Note that (u-u^) and (v-v^J define the pointing error 
The variance of the pointing error is then given by E(u-u^)^ and 
E(v-v^) . It should be noted that the mean pointina error becomes 
zero assuming that the phase errors added by the phasing tree have 
zero mean values. 


8.0 POWER TRANSFER EFFICIENCY 

Since the power efficiency is directly related to the beam 
forming and the beam pointing ability of the antenna, there are 
two different effects that degrade the power transmission efficiency, 
'hey can be generally classified under Random and Nonrandom effects. 
Figure 8.1 lists various sources of random and nonrandom effects while 
Figure 8.2 presents typical random effects added by the phase distribu- 
tion system. 

The combination of random effects has a mean value and a 
variance. The mean value of the random effects combined with the 
nonrandom effects produces the beam pointing error. In what follows. 
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Figure 8.1. Random and Nonrandom Effects Affecting the 
Power Transfer Efficiency. 
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we win Investigate the effect of the random and nonrandom parameters 
described above on the power transfer efficiency. 

8.1 $PS Power Transfer Efficiency Formulation 

8.1.1 Definition 

POWER TRANSFER Power Received by the 10 km Diameter Rectenna /„ 
EFFICIENCY * Total Power Radi ated~by the Spacetenna 


Figure 8.3 makes the idea clear. The power transfer efficiency can 
be redefined as 

POWER TRANSFER _ Power Output at Terminals A & B 
EFFICIENCY " Power Output at Terminals C & D 


This definition is convenient because the multiplying constants due to 
the propagation through the medium cancel out from the numerator and 
denominator. The computer program uses this efficiency definition. 

It may be noticed that the power received by the rectenna as well as the 
total power radiated by the spacetenna are functions of several 
parameters such as the spacetenna mechanical pointing error, jitters 
on the pointing error, jitters on locations of the radiating elements, 
feed current and phase jitters, to name a few. 

8.1.2 Computation of Efficiency 

Computation of total radiated power and the power received by 
the rectenna is done by integrating the averaged power pattern of 
the spacetenna over the appropriate regions. 

TOTAL RADIATED POWER = J’j' (averaged power pattern) sin ed 0 d 4 » 

H 
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where the averaged power pattern is obtained from the spacetenna 


configuration and the feed currents, and H is the hemisphere 
[e ■ (0 ,tt/ 2) and 41 » (0,2it)]. 


The actual radiating elements are the slots cut in the waveguides 


(slots are separated from each other by a distance of '/2, where \ is the 
wavelength of the waves in the waveguide). The pattern multiplication 


principle tells us that the radiation pattern of spacetenna is the multi- 


plication of the element pattern of the radiating elements and the array 


factor (based on the location of the radiating slots) of the antenna. 


Array factor of the antenna remains the same regardless of the nature 


of the radiating elements. The array factor of the spacetenna has a 


very highly peaked mainlobe having a 3 dB beamwidth of a few tenths of 


a minute of arc which is comparable to the angle subtended by the 


rectenna at the spacetenna center. The element factor of the slot 


is so flat around the boresight for several minutes of arc that 


the region of interest (the 3 dB beamwidth region) this factor 


could oe approximated by a constant. Hence the radiation pattern of 


the spacetenna in the region of interest, could be app*'Cximated by 


considering isotropically radiating elements instead of the slots. The 


spacetenna power pattern using the above approximation has two main 


parts comprising it. One part depends on the direction (f,,0 which 


will be called the anisotropic term and the other pan independent of 


(9,41) will be termed the isotropic part, i.e.. 


Averaged 

Power - E(f-f*) = Anisotropic Term + Isotropic Term 

Pattern 


The spacetenna averaged power pattern takes into account all the 
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system imperfections in the form of jitters. The isotropic term is a 
function of independent phase disturbances at the slot level. It should 
be noted that the isotropic tern carries into existance only if the phase 
disturbances at the slot level are statistically independent, e.o. the 
location jitter on the radiating slots this does not include, for 
example, the phase disturbance added by the phase distribution 
system because all the slots in the subarray are affected by the same 
phase disturbance. One more fact may be pointed out that if the slot 
element pattern is used in the computation, the only change is that 
the isotropic term becomes directional due to the element pattern 
and the total power due to this term reduces by the gain factor 
(of the slot) from the total power radiated by the isotropic term. 

Here we point out that 


THE TOTAL 
RADIATED POWER 



anisotropic term) sin eded(^ 


+ 



iostropic tern) sin edod()i 


(8.4) 


Figure 8,4 shows that the isotropic term is independent of (0,^), the 
direction (that is the reason it is called isotropic), lionce 


TOTAL 

RADIATED POWER 


^ (anisotropic term) sin Odediji 


Similarly, 

RECEIVED , 
POWER 


+ (isotropic term) Jj' sin 9ded(ji 


Constant 


JJ" (anisotropic term) sin eded(> 


e,iJ)eR 


(iostropic term) J j 

Q.^eR 


(8.5) 
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sin 9d6d^ 
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Components of the Power Pattern Due to Dipoles. 




where R Is the region of the rectenna. 

liumerical results pertaining to efficiency computation are 
given in Section 10.3. 

9.0 MODELING OF THE PHASE ERROR BUILD UP IN ’HE DISTRIBUTION OF PHASE 

STTHT TRrr 'STHummt 

The practical implementation of a retrofire array based on the 
phase conjugation principle requires the existence of a reference carrier 
at twice the frequency and constant phase at every phase conjugation 
circuit. The phase of this reference carrier must be exactly the 
same at every phase conjugation circuit; thus, the differential 
phase between the reference available at any two phase conjugation 
circuits is ideally zero. The distribution of this phase is a problem 
of considerable practical interest. The topic of this section is the 
modeling of errors in the dist*'ibution of this reference by the tree 
structure. 

9.1 Tree Structure for the Distribution of the Reference Phase 

As previously noted, the tree structure is a method of providino 
a constant reference phase to several physically remote points in 
space. This distribution is done in parallel; that is, several 
oscillators are "locked" simultaneously to one central reference 
oscillator. These remotely located oscillators serve as local 
reference oscillators to the oscillators located in their neighborhood. 
These new local references may now distribute the phase to another 
layer of "local" references. This process continues until it is 
possible to accurately distribute the phase to several phase 
conjugation circuits. This tree structure is illustrated in Fig. 9.1. 

(The single frequency method was drawn since itt simpler and the 
principle is the same.) A multiplied version of the reference drives 
four (in this case) PLL assemblies that are ultimately connected to 








the cables. The other end of each cable Is terminated in a short 
circuit or similar connection 50 that the signal transmitted down tne 

cable Is reflected back to the driven end of the cable. A portion 

of the Incident signal Is removed by the directional coupler and 

used as the reference for the next level In the tree. This circuit 

Is a remarkably straightforward solution to a complicated problem. 

It may be readily verified that the reference available at each side 

la-ld is exactly equal to cos(ut+e), the Input reference phase at level 

0 for an ideal implementation. Departures of this model from ideal are 

the topic of the remainder of this section. 

Keeping In mind that the purpose of the phase distribution network 
is to provide the same reference at all phase conjugation circuits then 
it is clear that the phase at references la through Id must all be 
the same. Any variations at this level cannot be removed by hardware 
at lower levels. First consider the transmission of the reference 
signal at the Input of the X2 multiplier at level 0 through to the 
reference regenerated at location Id. The signal is injected into 
the multiplier and then distributed by a power splitter (4-way 
for this example) and Injected into the PLL network. The PLL network 
generates a signal with the same frequency as the reference. The 
phase of this signal is equal to the reference phase minus the 
amount of phase shift that is contributed by the cable. Finally, 

the signal passes down the cable and is removed by the directional 
coupler. Since the hardware is designed to subtract out the phase 

effects in the cable only the hardware will contribute some phase 
shifts. These phase shifts are such that the references generated 
at level 1 outputs la-ld will differ from the reference phase injected 
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The points where phase shifts occur are at the multiplier, the 
four way power splitter, the PLL network (including the directional 
couplers, etc. that allow the two signals propagating in different 
directions to coexist on the cable), the reflective termination, and 
the directional coupler at the end of the cable that removes the 
desired signal reference at Id. For example, mismatches in the 
terminating Impedances in the power splitter will cause undesired 
phase shifts, and nonideal termination of the cable will result in an 
additional (undesired) phase shift of the signal returned to the PLL 
network. Each of these phase errors is a static error (neglecting 
thermal effects) and is set for all time. 

These static phase errors result in level 1 references with phases 
that differ from the reference injected at level 0. Since these phases 
are statistical in nature it is desirable to characterize thetn. Consider 
as an example the four way power splitter. Suppose there exists an 
ensemble of these four way power splitters. Each power splitter will 
have a specific phase shift between the input port and each output port. 
Given a specific output port then there is a certain mean value (averaged 
over the ensemble) of the phase shift between the input port and the 
chosen output port. In addition there is a variance associated with 
this mean. Now there is the possibility that the mean value (and 
associated variance) between the input and one of the other remaining 
output ports is some other value. This can be carried out for all 
the ports. The net result is summarized in Fig. 9.2. This is an 
extremely complicated model for a power splitter. The following basic 
symmetry assumption will now be used. If the power splUter is a 
symmetric design that there is no a priori reason to suspect that 
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Furthermore, the v« 'lances should al l be equa’ . This is 
probably a valid assumption since the Ideal power splitter has zero 
relative phase shifts between the output ports; thus any design that 
tends to favor a static relative phase shift between the output ports is 
a poor design. Thus the final model for a power splitter is to assume 
that the static phase shift between the Input and any particular 
output port Is a random variable drawn from some distribution with a 
certain mean and variance. This final model is illustrated in Fig. 9.3. 

This modeling philosophy can be carried through to the rest o^ 
the circuitry. For example now consider tne PLL networks, cable and 
termination. The mixers, phase shifters, power splitters and 
directional couplers will all contribute some phase errors. The net 
result is the phase of the recovered reference at a particular level 1 
node, say Id, will be different than the input phase at level 0. Now 
it nay be assumed that the relative phase shift betv;een the input to 
the PLL network and the recovered reference at a particular node is 
drawn from a particular distribution with a specific nonzero mean and 
specific variance. This means that the mean value dots not change 
from network to network! 

Then* are several final notes. First, it has been assumed that 
the hardware Is an ensemble with an associated ensemble mean and 
variance. This, of course, implies zero average static phase shift 
between the references recovered at level 1. There maj^ exist an a_ve_ra<^e 
phase shift between the reference at level 0 and aT[ the recovered 
references at level 1. Second, since any overall phase shift between 
the reference at level 0 and a1J_ the phase conjugation circuits will 
not squint the beam or produce any other deleterious effects, it may 
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as well be assumed to be zero. Thus, it may be assumed for the example 

considered here that the phase shift between the reference at level 0 

and any particular recovered reference at level 1 is a zero mean random 

2 

variable with a variance of o . Finally, note that the processes 
considered here are not ergodic. These phase shifts are static errors 
that will exist for all time. In spite of the fact that the ensemble 
average of the phase shift between level Oand level 1 is assumed to be 
zero, the actual phase shift will be some value that will exist forever, 
thus causing some static beam shift anomaly peculiar to that particular 


I arr^^y. 

9.2 Phase Error Build Up and Tree Correlation Matrix 

As we saw above, the phase distribution structure takes the 
form of a tree. Ideally the phase at any tip of the tree should be 
exactly the same as the master phase but as indicated above, there 
are sources of phase shift (random or otherwise) in which phase 
distribution structure which alter the value of the master phase 
delivered to the tip. Figure 9.4 depicts the sources of phase 
j shift per level of the phase distribution tree structure using 

MSRTS and the power splitters. The phase disturbance added at any 
poiiiC in the phase distribution tree travels through the tree towards 
I the tip giving rise to nonzero cross correlation of the phase at a 

tip with the phases at the other tips. The value of the cross- 
correlation depends upon the value of the variances of the phase 

errors added by the system to the signal as well as where in the 
distribution tree it was introduced. Figure 9.5 depicts a phase 
distribution tree having L levels of slave nodes, each node having 
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Figure 9.4. Sources of Phase Shifts Per Level of the Phase Distribution Tree 
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four branches. The phase error added by the phase distribution system 

at the 1 — level slave node comes from a distribution having a mean 

value and standard deviation a^. The subscripts on p and o indicate 

the north, south, east and west direction in which the slave node is, 

with respect to the center of the spacetenna. At the second level of 

slave nodes depending on which node is under consideration (north, 

south, east or west with respect to the immediate level 1 node) phase 

2 

error from the distribution having a mean value y and standard 
2 

deviation a is added to the phase error handed down by the immediate 

level 1 node. This combined phase error is then handed down to the 

third level slave nodes. The process of handing down the phase error 

from level to level continues until the tip of the tree is reached, 

at that point the phase with the accumulated phase error is used in 

the conjugation circuit to conjugate the pilot wave phase received 

at that point. An example of the phase error build up in a two 

level phase distribution tree having four branches per node is 

shown in Fig. 9.6. A general computer program was written for this particular 

application and it was executed with the value of phase errors and level 1 

and level 2 nodes from a Gaussian random number generator with zero 

mean and standard deviation of 5°. The resultant phase error at each 

tip is shown in the same figure. 

Thus we see that the correlation of phase at one tip of the phase 

distribution tree with the phases at the other tips depend on the 
routes taken by the signal from the center to arrive at those tips. 

An example of the correlation matrix for a two level phase distri- 
bution tree having four branches per node with the phase error 
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Figure 9.b. A Sample Phase Error Build Up in a Two Level 
Phase Distribution Network 







added at level 1 and at level 2 having a standard deviation of 5° 
(both the phase errors have zero means) is shown in Figure 
9.7. The matrix elements are normalized by ^ + f^Level 2‘ 

Note that the matrix has more zero entries than nonzero entries. An 
ideal matrix would be an identity matrix which results in maximum 
power transfer. 

9.3 Modeling of the Phase Error Build Up in the Distribution of 


^se ^^hen PAs are Present at Each Level of the Phase Distribution 
Tree 

Ideally the phase distribution tree transfers the master phase 
from level to level with no errors, hence it seems possible to put the 
PAs and the related conjugation circuits at each level reducing the 
cable weight necessary for the phase distribution system. When the 
circuit components are not ideal, randomness with nonzero mean (biases) 
are added to the phase being transmitted from level to level of the 
distribution system. Assuming that the master signal is a pure tone, 
the spectra of signals supplied to the PAs and the conjugation circuits 
is shown in Fig. 9.8. The spectrum seen by the conjugation circuit 
gets progressively broader and shifted as we proceed towards the tip 
of the phase distribution circuit. The variance of the phase after 

level is given by ^ ^Level i smaller than the variance 

of the phase at the tip of the phase distribution tree for total 

number of levels in the phase distribution tree. Thus the variance 
of the phase used in the conjugation circuits and PAs on the average 

becomes smaller than when the PAs and conjuoation circuits are at 

the tip of the phase distribution tree structure. This results in 
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higher gain of the spacetenna. The biases added by the phase distri- 
bution circuit components portray a different picture altogether. 

It Is well known that the conjugation scheme Is very sensitive to 
differential phase shifts of the currents feeding the radlaqlnt 
elements of the antenna. Fig. 9.8 shows that the differential 
phase shifts (biases) are considerable for the system with PAs and 
conjugation circuits at each level as compared to the phase distribu- 
tion system with PAs and conjugation circuits at the tips of the 
distribution tree having negligible differential phase shift. 

One serious source of biases is the power splitters. With 
some additional hardware at each power splitter the problem of 
biases can be eliminated. One such scheme is shown in Fig. 9.9. 

A VCO Is added to the power splitter assembly which has a feedback 
from the power splitter output, when the VCO locks. The biases are 
el iminated. 

10.0 REFERENCE FAR-FIELD SYSTEM PERFORMANCE AND TRADEOFFS FOUND 

vTAWARglH 

In this section of this report we shall demonstrate the far-fleld 
performance characteristics of the reference system phase 
control system as found using SOLARSIM discussed in Appendix IV. 

Since the problem is too massive and expensive to implement in the 
form of early experiments which could provide meaningful performance 
results, a computer simulation supported by analysis has been developed 
(Appendix IV). Effective use of simulation techniques as a system 
engineering tool allows one to establish performance tradeoffs as 
hardware technology develops without the extreme costs of full scale 
testing of various concepts. In tha following sections, certain 
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performance characteristics of the SPS reference system are presented 
as found using SolarSIm— an analytical simulation of the SPS phase 
control system. 

10.1 Phase-Error Build Up and Tree Structure Optimization 


In order to evaluate different phasing tree structure, one has 
to study the effects of the tree structure on the spacetenna perform- 
ance. In particular, the boresight gain reduction and the spacetenna 
pointing error are two Important performance measures for merit 
comparisons. The boresight gain reduction Is proportional to the 
transmission efficiency of the microwave transmission system. The 
pointing error measures the ability of tie spacetenna to focus on the 
center of the rectenna. SOLARSIM has been developed to In'^estigate 
these effects. The current reference phase distribution tree consists 
of four levels and Is shown in Fig. 10.1. Note that a total of 102400 
nodes (PAs) at the last level are fed by the structure. Phase noises are 
Introduced to the nodes at each level by the MSRTS discussed in Section 
4.2.1. In addition, In the last level (PA level), additional phase 
noises are Introduced by the power transponders . 

10.1.1 Effect of Tree Structure on Main Sean. Gain Reduction 

The main beam gain reduction Is a function of the magnitude of 
and the level where the phase errors are introduced by rhe phasing 
tree. It Is interesting to see the extent to which the phase errors 
introduceo at each level Individually contributes to roduce this g^ln. 

In Fig. 10.2, we show the effect of ♦’he phase errors introduced at a 
particular level In the phase distribution tree, while the phase error 
at the other levels are neglected. A phase jitter of 5” is shown for 
a four level tree. The normal main beam gain loss of gain in dB s 
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Figure 10.’, Four Level Phase Oistribution Tree. 



Is plotted In Fig. 10.2. 

The quantity is obtained when all the phase errors Introduced 
in the phasing tree are set to zero. From Fig. 10.2, we see that as 
far as the gain reduction is concerned, phase errors introduced in 
the last level (4^^ level) results in the worst performance. Intuitively, 
phase errors introduced at the beginning levels of the tree cause 
the total phase error build up at the radiating elements to be 
correlated. As a result, individual patterns all add up in the same 
general direction. On the contrary, when phase errors are introduced in the 
last level, radiating elements have uncorrelated phase errors which 
results in reduced main beam gain. In any case, the differences are 
small . 

10.1.2 Effect of Tree Structure on Rf»1S Pointing Error 

The ideal spacetenna attains its peak power density in the direction 
of the rectenna center. Under the assumption that the power beam is 
focused properly, the purpose of the phase control system is to phase 
its individual radiating elements so that the boresinht gain is 
maximized. When random errors are introduced by the phasing system, 
the antenna power pattern peaks at some other direction other than 
the boresight. The variance of this direction is spacetenna pointing 
error. 

Depending on the level where phase errors are introduced, the 
resulting pointing error is different. From Fig. 10.3, we can see 
that the phase error introduced at the first level causes the largest 
rms pointing error while the phase error from the last level (level four) 
has the least impact. The same trend is observed from the 8- and 9-level 
phase distribution tree which is depicted in Fig. 10.4. In essence, if 





Figure 10.4. RMS Pointing Error as a Function of the Level where 
Phase Errors are Introduced. 
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Indtpendent phase errors are Introduced at the last level, the resultant 

phase error variance on the beam Is reduced roughly by the number of 

radiating elements (101552). This follows from the law of large 

numbers. However, If the phase errors are Introduced at the first 

level, the resultant phase error variance on the beam Is reduced by the 

ntmiber of nodes In the first level. For the baseline (4-1eve1 tree), 

that number is 25. The square root of the ratio 101552/25 1s roughly 

18 dB, which can be observed In Fig. 10.3. 

10.2 Performance Evaluation Using SQLARSIH 

10.2.1 Phase Error Build Up Effects on Far-Fleld Power Vattern 

One of the current capabilities of SOLARSIM Is to compute and 

display the received power pattern. The set of patterns shown in 

Figs. 10.5 to 10.9 represents Its use to investigate the effect of 

phase error on the received power pattern for the 4-1eve1 distribution 

tree baseline. Fig. 10.5 is a typical power pattern with the phase 

jitter per level set to be 5®. The power pattern is normalized so that 

it is unity if all phase errors are negligible. This corresponds to a 

2 

maximum microwave power density of 23 nW/tm at the center of the 
rectenna. The power pattern is shown as a function of e which is the 
angle from the boresight of the spacetenna. At any point on the rectenna, 
its distance d from the center is approximately given by d = R6 where 
R is the distance of the spacetenna from the rectenna (3.7x10^ Km). 

The edge of the rectenna is approximately equivalent to half a minute 
of arc, which is depicted in the figure. Both the power patterns for 
the 0 dB taper (dashed line) and 10 dB taper (solid line) are shown 
in Fig. 10.5. We can observe that the power taper serves to bring down 
the sidelobes while broadening the main beam. For the reference system 



m 


- 95 - 











4 LEVEL PHASE DISTRIBUTION TREE 
SIGMA PER LEVEL - 180* 


1 0 dB TAPER 
... 0 dB TAPER 
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(10 dB tiptr). the pOMtr dtnsity 1i 6am to -13.3 dB at tht rtcttnna 
td 9 «. This eorraspofids to a poMtr dansity of 1.08 mW/cn^. Tha ttrlts 
of figurti fhoM that as tha phasa jittar par laval Incraasas, tlia 
borasight powar dansity dacraasas and tha sidaloba lavals bacoma 
closar togathar. Finally, for an rms phasa arror of 180* par laval, 
l.a.. Initially no phasa control, tha spacatanna appaars to ba an 
Isotropic radiating alanant (saa Fig. 10.9). Undar this condition, 
tha po»i«r dansity Is approxlmataly 3 x 10"® wM/cw^. This laval Is 
far balOM althar tha U.S. radiation limit (10 mW/cm^) or tha USSR 
guldallna (10*^ mU/an ) and hance does not appaar to posa any 
biological hazard. 

10.2.2 Main Baam Sain Raductlon Qua to Phasa Error Build Up In Trae 
0? Phase Control Systa>:> 

Wa have seen qualitatively how the efficiency of the microwave 
transmlsslm system Is reduced If the phase error build up 1$ not kept 
under control. Fig. 10.10 Investigates this effect by plotting the 
borasight gain reduction vs the total rms phase error build up through 
the phasing tree. As we have seen earlier, the shape of the power 
pattern 1$ relatively Insensitive to the phase error, at least to the 
first few lobes shown provided the rms phase error Is not excessive 
(e.g. < 30*). In that case, the main beam reduwtion Is proportional 
to the power received by the rectenna. In Fig. 10.11 and In what 
follows, total rms phase error Is tefined to be the total phase 
error build up through the phase distribution tree, assuming each 
level In the tree adds In the same amount of phase error. The phase 
error processes Introduced at different nodes on the same level are 
assumed to be Independent of one another. For an efficiency of 96%, 

Fig. 10.11 Indicates that one can only tolerate about 11.8* of total rms 
































phase error. The allowable phase error build up In the tree structure 
for a fixed efficiency Is tabulated In Fig. 10.11. For the 96 % 
efficiency, the total (RF) phase error build-up of 11.8® reflects into an IF 
phase error build-up of 2,4® for the 490 MHz IF baseline. For a four 
level tree (baseline system) the maximum error tolerable Is 5.9® per 
level. If a nine level tree Is used, we have to tighten the phase 
error specification to 3.9® per level. In this respect, the 4 level 
tree system seems to be superior. However, since different amounts 
of hardware and layout Is required for different structures, the above 
remark Is not conclusive. 

As we have eeen earlier, the phase errors Introduced at the last 
level result In the most severe reduction In the main beam gain. This 
effect Is further Investigated In Fig. 10,12 in which the total phase 
error build up through the phasing tree up to, but excluding the last 
level. Is held fixed at 15® rms, and we vary the phase error Introduced 
at the last level. The resultant gain reduction is Identical for the 
four and the nine level tree. In general, one expects the magnitude 
of the phase error Introduced In the last level to be higher than 
the earlier levels, since we must now take into account the phase 
errors introduced by the power transponder, in addition to the phase 
reference distribution system. 

10.2.3 RMS Pointing Error Due to Phase Error Build Up in Tree of 
Phase Control System 

The sensitivity of the rms pointing error to the total phase error 
build up is investigated in Fig. 10.13. As to be expected, the pointing 
error is directly proportional to the phase error. It is of interest 
to relate the rms pointing error to the probability that the beam misses 
the rectenna center by a given amount. This is rather straightforward 
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If Me assume that the pointing error Is approximately Gausslanly 
distributed. For a typical total rms phase error of 15®, we can see 
from Fig. 10.13 that the rms pointing error Is 0.019 min. This 
translates Into a probability of 0.99 that the main beam will be 
restricted to 0.03 min off boresight. In terms of distances, this 
means that the probability that the main beam misses the rectenna center 
by 300 meters Is less than 10 . 

10.2.4 Sumnary of Performance of the Baseline Phase Control System In 


Terms ofMain Beam Gam and Pointing Error 


The table In Fig. 10.14 summarizes the performance of the baseline 

phase control system. The projected efficiency Is compiled as a function 

of the RF phase error build up (both total and per level). The associated 

rms pointing error and the 99 percentile range of the main beam are also 

given under the corresponding phase error conditions. For a target 

efficiency of 96%, the RF rms phase error allowed 1s 12® total or 4® 

_2 

per level. The corresponding rms pointing error is 10 min while the 
99 percentile distance is 437 meters. Hence if we can keep the phase 
error to be less than 4® per level, the performance of the phase control 


system Is satisfactory. 

From [11], the gain reduction is approximately 

2 2 

exp(-o^) where a. is the total phase error variance and the rms 

9 9 

pointing error is proportional to the rms phase error. These 
relationships were developed for a linear array. Nonetheless, they 
also seem to apply for the present system. As a matter of fact we 
can verify that G/Gg » exp(-a^) and * (?^/4.5xl0^ from Fig. 10.14. 
10.2.5 Effects of Antenna Imperfections on SPS Performance 


In an ideal retrodi recti ve array where phase conjugation is performed 


cJHinCoi 
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on each and every radiating element* the performance of the array Is 
Insensitive to Imperfections on the array structure. In reality, the 
phase conjugation required for retrodirectivlty Is feasible only at 
the subarray level. As a result, the performance of the spacetenna 
as a retrodi recti ve array Is severely affected by the layout, size, 
surface roughness and tilt of the subarray. The latter phenomena 
are caused by structural Imperfections due to engineering tolerances, 
bending and shearing forces, aging, and thermal effects, to name a 
few. SOLARSIM provides the capabilities to Investigate these effects 
on the spacetenna performance. 

10.2.5.1 Spacetenna Layout, Tilts and Location Filters 

The 10 dB power taper In the present reference system is achieved by 
dividing the spacetenna Into 10 distinct power density rings. A 
quadrant of the spacetenna Is shown In Fig. 10.15. The radii of these 
rings run from about 94 meters to a maximum of 500 meters as shown. The 
power density per ring Is varied by changing the number of PAs feeding 
per unit radiating surface. For example, the number of PAs feeding a 
10.4 m X 10.4 m square of slotted waveguide changes from 36 to 4 from 
the inner most to the outer most ring. The power taper from the 
spacetenna center to Its edge Is thus 9.5 dB co be exact. 

The current reference system assumes that phase conjugation 1s performed 
on each PA that feeds a subarray. By a subarray we mean a group of 
slotted waveguides grouped Into a square In the center of which phase 
conjugation Is performed. Hence only the central radiating element Is 
truly retrodi recti ve. The other elements on the same subarray will 
experience differential phase shifts defined by the projections of the 
location of the elements along the boresight. The size of subarray in 
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9 2 

the reference system varies from 1.73x1.73 m to 5. 2x5. an from the center 

to the edge of the spacetenna. The details are given in Fig. 10.16. 

Also shotfKn in the same figure is the fixed size subarray geometry 

conceived by Boeing. In that configuration all subarrays have a 

2 

fixed size of 10.4 x 10.4 m . Only one conjugation center is provided 
on each subarray. 

For a perfectly flat spacetenna, the normal directions of each 
subarray are all identical, i.e., they all coincide with the normal 
direction of the spacetenna as a whole. Since the spacetenna is never 
perfectly flat in reality, each subarray will point in a different 
direction. For SOLARSIM, we adopt the subarray tilt model depicted in 
Fig. 10.17. The dashed line shows the ideal position of the subarray. 

With tilts, the subarray is depicted by the solid line. Note that the 
tilted subarray makes an angle e with the x axis and an angle 6 with 

X y 

the y axis. The angle between the present normal direction and the ideal 

2 2 2 

direction is given by e ■ The tilted angles e and e have mean 

y X X y 

components denoted by E(e ) and E{e ). They represent the errors in the 

X y 

SPS steering mechanism made in pointing the spacetenna towards the 

rectenna center. The mean values E{e ) and E(e ) are the same for all 

^ y 

the subarrays and they are referred to as the mechanical pointing error 

in the x- and y-direction. The perturbation about these means, i.e., 

0y-E(0y) and 0„-E(e„) represent the variations of the tilt angles from 
XX y y 

one subarray to another. They are assumed to be Gaussian with variances 

2 2 
0 and 0 . 

X y 

In the above model, we have modeled the subarrays to be flat surfaces. 
This assumption is Justifiable since we are only interested in the global 
behavior of the subarray. However, when we are concerned with the 
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Figure 10.16. Variable and Fixed Size Subarray Gecxnetry 
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Figure 10.17. Effect of Subarray Mean Tilt and Jitter. 



radiating elements, we need a more refined model for SOLARSIM. Figure 


10.18 defines the location uncertainty of the radiating elements and 
the conjugation point primarily resulting from the unevenness of the 
subarray surface. The average plane through the subarray defined the 
Ideal subarray used to model the tilt effects. However, the actual 
locations of the radiating elements and the conjugation point will be 
perturbed from the average plane. A cross section of a typical subarray 
surface running through the conjugation point parallel to the y>ax1s 
Is shown In Fig. 10.18. The location jitter Is defined to be the 
projection of the location of a point along the line joining the Ideal 
center of the subarray and the center of the rectenna. 

10.2.6 Sensitivity to Tilt Effects as a Function of Conjugation Level 
In the absence of tilt effects, the performance of the spacetenna 
Is not sensitive to the conjugation levels. In practice, the presence 
of tilt effects determines the level required. Fig. 10.19 compared the 
sensitivity of the gain reduction G/6 q vs the mechanical pointing error 
'With the location jitters (assumed equal for conjugation points and 

n — T" 

radiating elements) as a parameter. A typical rms tilt o. ■ /a* +o 

V Ojj By 

of 2 min Is assumed. We can observe that the baseline system degrades 
more gracefully than the fixed subarray (Boeing) scheme. 

10.3 SPS Power Transfer Efficiency 


In Section 8 the power transfer efficiency was defined. That definition 
Is used to generate the results below. 

The computer program evaluates the received power as given In 
the equation (8.6), I.e., ♦c(0,2ir) and 0c(O,O.464' ). The total radiated 
power Is computed In two parts. Tho first part Integrates the anisotropic 
term of the power pattern over the main lobe and 5 sidelobes, I.e., 'H* 
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Figure 10.18. Jitters Added Due to Location Uncertainty of 
Radiating and Conjugation Points. 
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In tquatlon (8.5) Is approxlnattd by ac(0,2w) and ec(0.ma1n loba * 

5 sidalobts). Tht stcond part computts tha Isotropic part of tha 
poMar pattam and multlpllas It with tha ‘constant* sho'/m In aquation 
(8.5). Notica that tha constant Is hand computed with H to ba 
ac(0,2ir) and ec(0|W/2). Addition of these two parts gives the total 
radiated power. The In^ortant thing to notica Is 

. -2w j/2 2iu 

(Isotropic term) / / sin ededa » (Isotropic term) j j sin edsd# 

•'a-0*'e-0 a-*0 6c(5B.L.) 
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* -t/2 

I sin ededa » 


a*o 0 "O 
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sin edoda 


a-0 ec(5B.L.) 


(10.3.1) 


(10.3.2) 


The computer program computes the received power and the total radiated 
power as described above and their ratio gives the efficiency. This 
efficiency Is plotted against the total rms phase error In F1gurel0.20. 
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10.4 Observations and Improvanent s 

The first observation one makes from Fig. 10.20 Is that the 
efficiency of the system Is below 952 even with no phase errors present 
as Indicated by curve ^ . At a phase error of 10® rms, the efficiency 
drops to about 902 with no other error present. The following suggestions 
are made to increase the efficiency: 

a. Put more conjugators in each subarray (increase in 
efficiency would be minimal). 

b. Use radiators with more directivity instead of dipoles 
(slots). 


c. Increase the size of the rectenna. 
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ANALYSIS OF THE TWO-FREQUENCY PILOT SIGNAL CONCEPT 
AND AMBIGUITY RESOLUTION 


As we saw above, the necessary evil of frequency separation in 
the single frequency separation in the single frequency conjugation 
scheme introduces the phase correction necessary and that necessitates 
the computer. Besides there is a question of feasibility and 
accuracy of such a system. There exists, however, an easier way to 
circumvent the frequency separation problens. Consider an array of 
transponders located at positions Rq,Ri ... . ,R j^ relative to the location of 
the pilot transmitter. A carrier Uj(t) » e^“^ emitted by the 
pilot is received at the i— transponder. This received signal 
is The amplitude factor has been neglected 

for the sake of this illustration. The net phase shift 
accumulated by the signal in traversing the distance |R^| is 
2irjR^|/x. The complex conjugate circuit conjugates the amplitude 
factor and the reradiated signal is 

This reradiated signal is received at the pilot location. This 
signal is 2nl^il/x), gja»t, This is independent 

of the location R^. Thus the signal received at the pilot tone 
location is the same for any transponder location. Since the 
reradiated signals from each transponder are all in phase at 
the pilot transmitter location then the power is maximized 
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and the peak of the array beam Is In exactly the same direction 

f’ 

I as the pilot tone transmitter location. 

I insider now the case where the down link frequency Is 

I Then the reradiated signal Is and the 

received signal Is 2 «|^^|/Xq)^ The phase Is 

1(1/X-1/Xq). Since ui^uq, Xj'Xq and now the phase does 
depend on }R^|. Now the beam does not peak In the direction 
^ of the pilot signal. Note that -2ir|it^ |(l/x-l/Xg) * -2wli^^ |.( w-uq)/c 

The amount of the phase shift Is proportional to the difference in 
frequency between the uplink and the downlink. This suggests that 
if two frequencies were used for the uplink, o»q* wq+Lu and ug* wq-Au, 
then the phase shift for each frequency could be measured and the two 
phase shifts could be summed and divided by two. This average 
phase shift Is just the phase shift that will be undergone on 
the downlink carrier with frequency Ug. 

Phase shift at u » ^ ug+A« • ~ |R^| (wg+au) 

Phase shift at « » U 2 ^ “O"^ " ^ 1*^1 1 
2»|RJ 

Average Is [(ug+au)+(ug-au)]/2 

2irlRj 

• ~T“ 1) 

* phase shift on downlink 
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The phase shift on the two frequencies Is measured and averaged. This phase 
shift Is added to a carrier at the center frequency wq*(ij^ - tu^/ 2 At this point 
the remaining operation Is the same as an Ideal Van Atta array, 

I.e., the complex conjugate of the signal formed and this 
signal Is reradiated. This operation Is Illustrated In Figure 1.1 
The uplink signal Is 


r(t) ■ d(t)[cos{a»^t+0^)+cos(u2t+82)] 


(I.l) 


where u^*wq+&u, « -^1R,-1(uq+Aw), * c ^^iI 

and d(t) Is some data. 

eo“(0i+02)/2 cosjuQt) 

Phase Phase cos((ijQt+0Q) 


Phase 


Phase 

Estimate 


Modulator 



Complex Conjugate 
^ cos 2u„t 




C0S(uQt-6Q) 


Figure 1.1. Block Diagram of Two Frequency Phase Estimation 
and Conjugation Transponder. 


The problem of performing the phase estimation and phase 
modulation may be solved fairly easily. A circuit that 
performs this Is Illustrated In Figure 1.2. 








Loop Filter 


sin 2(uQt+0Q) 


+sin(uQt+eQ) ^ vC/ 

Figure 1.2. Phase Estimator and Modulator Circuit. 

This circuit operates as follows: 

2 

x(t) ■ r^( ) • Ccos(u^t+0^)+cos(u2t+92)3 

■ cos^(u'jt+0^ )+cos^(<jj2t+02)'^2cos(ii»it+e^ )cos(u2t+02) 

( 1 . 2 ) 

The desired signal is contained in the product terms 

2cos(u^t+0^)cos(u2t+02) • cos 2{awt+A0)+cos 2((DQt+0Q) (1.3) 

»,here A0 - % ande^-^ll^^i wq* The second term is 

tracked by the phase locked loop. The PLL output is sin 2(uQt+0Q). 
This is frequency divided by two in a frequency divider. The output 
is +sin(ugt+0Q). Except for the +(0*-180“) phase ambiguity due 
to the divide by two circuit this signal is exactly the desired 
signal. It has the correct center frequency Uq = (u^+u2)/2 and 








the same phase shift 6 q ■ 2v|Rf|uQ/c that vMuild have occurred 

If the uplink signal was a single sinusoid at frequency uq. 

One serious problem remains to be solved. The divide by two 

results In a 0*-180* phase ambiguity. This Is a potentially 

• • 
disastrous situation since phase reversals represent 

cancellation of signal power and thus result In a serious 

deterioration In array performance. This phase ambiguity may 

be resolved by the following circuit which Is connected to points 

A* B and C In Figure 1.2. 


PLL 


1 

• 2 

cos 2{&ut+AB) 




+COS(u„t+en) 



+d(t)cos(awt+Ae) 
Figure 1.3. 


The circuit works In the following fashion. At point C Is available 

« 

the reference ls1n(<DQt+0Q) will contain the 180® phase ambiguity. This 
is delayed by X^Z 90® to get +cos(u)Qt+eQ). This Is 
multiplied by the input from point (a). The result of this 
multiplication is jfd(t)cos(Aut-^Ae). This contains the data on 
a carrier at the difference frequency. This signal contains the 
180® phase ambiguity. At point ^ is the second harmonic of 
this carrier which Is the first term In equation (1.3). This 


h. 
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coi|>onent Is removed using the PU. The output of this PLL is 
frequency divided by two to give ^cos(0(iit-«-Ae) where the * 

Indicates another 180* phase ambiguity due to this second 
divider. These ambiguities may now be resolved using the 
data d(£). This Is a two step acquisition procedure. First 
d(t) ■ +1, for a fixed period of time. During this time the 
output +d(t)cos(^t’i'Ae) is compared with ^os(out>ae) at the 
output of the divide by two. These two signals are compared to 
determine whether or not they are in phase or out of phase. 

This Is used to resolve the ambiguity in the circuit of Figure L3. 
The phase compare circuit is relatively straightforward and 
is illustrated in Fig. 1.4 below. Thus the output of the 

;Z Output 


0 Out of Phase Out of Phase -♦Advance 



■ (2A) In Phase In Phase -*>Do Nothing 


Figure 1.4. 

PLL/ 12 combination can be phase correctly with the data 
bearing signal. This signal is now used to denwdulate d(t) 

as illustrated in Fig. 1.3. The data is now modulated with a sequence 

of I'-s and -I's. If the number of Ts is different than the 

number of -I’s over a given interval of time then the number 

of I's may be counted and used to resolve the final phase ambiguity 

in the divider in Figure 1.2. 





Th# squaring typ# of circuit In Fig, 1.2 contains a divider 
that Is probably Impossible to Implement at microwave frequencies. 
A Costas loop may be used Instead which eliminates both the 
squarer and the divider. This removes one problem but 
Introduces the problem of differential phase shift In the two 
arms of the Costas loop. Either one of these circuits must 
operate on the entire range of frequencies of Interest, I.e., 

This Includes the center frequency and thus the 
possibility exists for leakage of Uq terms from the antenna back 
through the diplexer to the receiver. This Isolation question 
Is crucial to the operation of these circuits. notch may 
need to be Inserted at the center frequency. 



APPENDIX II 

SURVIVABLE TOPOLOGICAL CONFIGURATIONS FOR THE 
SPS PHASE CONTROL SYSTEM 

Introduction 

Th« object of this appendix Is to apply the mathematical theory of 
graphs and topology to the problem of evaluating the reliability 
and survivability of the SPS Phase Control System due to failures. 

Graceful degradation in the power transmitted is required as failures 
occur In the Phase Control System. In doing so the cable lengths, 
number of connectors, etc. will be minimized for an a priori 
specified connectivity resulting from failures of parts In the Phase 
Control Centers and/or the Power Amplifiers and Interconnecting cables. 

It Is one purpose of this report to develop optimum survivable topological 
configurations which conform to a set of prescribed performance and 
structural measures. It Is desirable to derive PCC and PA configura- 
tions which require the least overall cable length while attaining 
a properly low value of m ■ minimum number of cable pieces and connectors. 

The synthesized PCC and PA configurations are required to be 
survivable to a specified degree when failures occur. These configura- 
tions should perform effectively under failures of certain PCCs and/or 
PAs or Interconnecting cables. We will measure the survivability of 
the Phase Control System In terms of Its Invulnerability to PCC and/or 
PA failures, I.e., we required grac»ful degradation of power trans- 
mitted as a function of failures In the Phase Control System. We 
note that the resulting survivable PCCs and/or PAs are always Invulner- 
able to cable failures. A variety of PCC/PA configuration Invulnerability 
measures will be used Including deterministic and probabil istic invulner- 
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ability measures. TtM results produced Include the ability to predict 
PCC and PA Interconnectivity so as to optimize power output as a 
function of failures In the SPS Phase Control System. 

In what follows a more detailed discussion of the above concepts Isqlven. 

This Is followed by a few examples with recommendations for further Investloatlon 
II. 1 Topological Structures and Performance Configurations 

The SPS reference phase distribution systems can be modeled 
from a topological configuration point of view, as follows. 

Each amplifier node (containing the associated radiation antenna 
and other processing elements) Is taken to be a node (vertex) 

In an underlying graph. We thus Identify (n>1) nodes, 

Vp ■ {v^ .Vg,... ,v^} . 


We also Include In the graph a central node, denoted or y^. 

■ The central node v^ distributes phase Information to the "regular" 
nodes called Phase Control Centers (Refs. 1,2). Each regular node must 
be able to obtain this "reference" phase Information for It to effectively 
survive. A similar situation will exist for t^e power amplifiers. 

To distribute phase Information from the central (master) node to 
all others, one must connect by "cables” the regular and central nodes. 

A cable connection between node v^ and node Vj Is designated by a line, 

denoted (v^.Vj) or (1,j), or e^j. The collection of nodes V » {v^,1»l n} 

and of lines E connecting them, 

E * 1e^j,1,j”l,...,n }, 

form the underlying graph G ■ (V,E). This graph describes the pure 
topological structure of the network constituting the PCCs and the 
cables Interconnecting them. 
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To Incorporate geographical distances Into ;i^e typological 
description of the antenna network, we generate the weighted graph 

««• 

Here we have described the line connections by a distance matrix W. 
where 


W ■ {«^j, 1,J“l,2,...,n} 
and 

■ distance from vi to Vj, If they are connected 

■ •, otherwise. 

Thus, W describes the network line connections as well as the associated 
distances of the connections. 

It Is the purpose of our study to develop optimal survivable 
topological configurations Gy, which conform to a set of prescribed 
performance and structural measures. The latter constitute a proper 
set of constraints to be Incorporated In the tree synthesis 
procedure. We outline In the following, the najor such topological 
constraints. 

. 1. Number of Lines and Overall Distance Measures 

The overall number of lines in the underlying graph is denoted by m. 
Thus, m represents the number of cable pieces used. Furthermore, since 
each cable piece Is connected at Its ends to two nodes, then 2m is 
equal to the number of connections (and connectors) used In the system. 

It Is therefore of Interest to develop a survivable topological 
configuration which will have the minimal number of lines (minimum m) 
and thus the minimal number of cable pieces and connectors (as well as 
local connections). 
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If the cable pieces to be used In the network are, generally of 
equal (or nearly equal) length, the number of lines m also represent 
the overall length of cables used In the system. Mini nizinq m, would 
thus also entail the utilization of the minimum overall cable length. 

If cable pieces are of varying actual lengths, the overall cable 
length Is given by 

i“i ■ EE “ij 
i J 

It Is then desirable to derive a network configuration which requires 
the least overall cable length |W|, while a;ta1n1ng a properly low 
value of m • minimum nun^ber of cable pieces and connectors. 

1 1 . 1 . 2 . Survivability, I nvul nerability and Connectivity Measures 

The synthesized tree network Is required to be survivable. The tree network 
should perform effectively under the failures of certain regular nodes 
or lines. Node failure Impose more serious effects on the network 
performance. They are also more liable to occur (than llne-cable- 
fallures) In our situation. Me will thus measure the survivability of 
the network configuration In terms of Its invulnerability to node 
failures. We will note that the resulting survivable netw^ ks are 
always also Invulnerable to line failures. 

A variety of tree network Invulnerability measures can be utilized 
They can be classified Into two general classes: deterministic and 
probabilistic Invulnerability measures. We will mention here a set 
of deterministic and probabilistic connectivity measures which are 
proper for the tree network. 
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Network Connectivity Measure . 

A netwcfk of PCCs or PAs Is said to be k-connected if and only if 
the network remains connected under the failure of (k*l) or less 
nodes; k=l , 2 ,... . 

A network is said to be connected if and only if there exists an 
least one path betwee- any two nodes; thus, any two nodes can "communicate" 
to one another phase nifcn:.?tion. 

It can be shown that if a network is k-connected then there exisU 
at least k (vertex disjoint) paths between any pair of nodes v. and v.. 

* J 

Thus, there exist at least k paths between v. and which have in 
common only nodes v.,v^. Then, upon the failure of any node and its 
associated path, the remaining paths could be used for intercommunication 
between the corresponding terminal nodes. 

For our phase distribution network, a k-connected topological configuration 
ensures at least k (vertex disjoint) paths between the central node 
and any regular node. Subsequently, the central node can conmunicate 
its phase information to any regular node even if (k-1) or less other 
regular nodes have failed. 

Note that a 1 -connected network can assume a tree structure (having 
no cycles). Upon the failure of any single regular nodes, certain nodes 
become disconnected from the center of the tree. 

A 2-connected network can be realized by a cycle structure. 

We will show how k-connected networks with minimum number of lines, 
any k>l , satisfying certain prescribed topological constraints can be 
developed. 


j 


J 


JUinCo m 


mk 


J 


- 133 - 



I I ? 
i 

f 


I ! 

f 


■JlunCc 


om 


A Probabilistic Connectivity Measure . 

We can regard node failures as a random phenomena. The latter is 
then statistically characterized by the joint distribution of failure 
events associated with the network nodes. 

For example, one can assume node 1 failure to be statistically 
Independent of node j failure, 1?*j, with a given probability that 
node 1 will fall within a certain period of time. 

A proper probabilistic invulnerability measure would then be given 
by p, where: 

p = min Pr{node v. Is connected to the center} . 

1 : 1>2 

Thus, If a network G has p(G) * p, then the probability that a reqular 
node will be connected to the center Is not smaller than p; It will be 
disconnected from the central node with a probability not higher than 1-p. 

We will show how measure p can be computed and incorporated into 
the topological design. 

11.1,3. Nodal Degree Constraints 

The degree of a node is defined as the number of lines (cables) 
connected to this node. 

Practical considerations Impose actual constraints on the degree 
of a regular node, as well as on the (usually higher) degree of the 
central node. In our development of optimal topological configurations, 
we will Incorporate these Important degree constraints. 

We note that If a network is k-connected, the degree of each node 
must be at least equal to k; i.e., 

deg(v^) ^ k, each i. 
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It Is generally desirable to synthesize the antenna topology 
such that every regular node will have a degree not much higher than k, 

k < deg(v^) < k + N|^, each i ^ 2, some N|^ . 

The central node should have a degree not higher than some constraint 

bi b > k, 

k <. deg(v^) £ b . 

Note that if all network nodes have a degree equal to i£k, then 
the network has a number of lines m equal to 

m ■ ji-n, 

where n is the number of nodes in the network. Thus, increasing nodal 
degrees, clearly induced an increase in the network number of lines m. 
Also, increasing the network connectivity k, implies higher values for 
nodal degree, and a higher number of network lines m, as 

1 

m = 2^ deg(v^. ) . 

i=l 

Such considerations will be an integral part in our topological studies. 
.4. Network, Radius and Radius-Stabi ’ ity Constraints 





For our antenna network, we are Interested 1n defining the 
following radius r ■ r(G) measure: 

r » r(G) » max d(v.,v,) . 

i: i>2 ^ ' 

Thus, for our centralized network G, the radius r(G) is equal to the 
maximal distance of any regular node from the central node. 

The phase information distributed from the central node to the 
regular nodes across the network cables, experiences phase delays along 
each cable. To limit the overall phaFt ielay experienced by the phase 
stream, we need to limit the lengths of the paths used for its distri- 
bution. Thus, we impose a constraint on the maximal radius of the 
topological structure: 

r(G) _< rQ, 

for same proper rp value. Subsequently, any phase message between the 
central node and any regular node, is routed through a path not longer 
than rg cables (lines). 

One can also incorporate actual geographical distances to define 
the above mentioned radius, diameter and distance measures. 

Furthermore, from the viewpoint of graceful degradations of 
efficiency under failures, it is important to derive a topological 
structure which is radius-stable . For such a configuration, upon the 
failure of (k-1) nodes, or less, the radius of the antenna remains 
unchanged, or just changes by not more than a prescribed value. Thus, 
upon failures of certain nodes in the antenna network, we not only 
provide that it remains connected, but also prevent the new distances 
between any regular node and the center from increasing arbitrarily. 
Subsequently, as nodes fail, the remaining nodes are still able to 


i 
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receive phase Information which does not contain much extra phase 
delays. 

We note that In making the PCC/PAs k-connected we have designed 
a topological structure which sets (at least) k (vertex disjoint) paths 
between v^ and any regular node. To make the network k-connected and 
radius-stable, we will also constrain the latter k paths to be not (much) 
longer than the prescribed network radius Vq. In this way, we will be 
able to provide acceptable phase distortions even under nodal failures. 

We have recently developed new procedures and graph theoretical 
techniques, for the synthesis of optimal survivable topological structures 
applicable to PCC and PA layouts. We will extend these results and 
techniques to apply to the present problem and thus obtain the desirable 
optimal survivable topological configurations for the SPS phase control 
system. 

In the following we Illustrate the underlying topological design 
problem by presenting various optimal topologies of interest. As noted 
these structures will be extended and applied to satisfy other more 
general constraints that might be incorporated into the antenna network 
design problem. 

II. 2 The Reference Phase Distribution Network and the Graph Adjacency 

In the performance analysis of the reference phase distribution network, 
the following connectivity matrix has been defined: 


F - {^ij*^»J “ U2,..., n} 

where is a submatrix (row vector) of order Ixn, defined as 


F^j * {0,0,..., 0} if Ifj, each i,j , 


F.^ » ^*iT®i2’-” ’*in^ ’ ’ ‘ 
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where, for 1j<j: 


and 



1, if vertex v. Is connected by a line to 
vertex Vj ’ 

0, otherwise 



Thus, row vector represents the connections made between node v^. 
and all other nodes. Matrix F subsequently describes the topological 
structure of the cable network interconnecting the power amplifiers 
(nodes). 

Clearly, for the sole sake of the network topological description, 
matrix F can be reduced to the adjacency matrix of the underlying 
graph . 

For a graph G = (V,E), with set of nodes V = {v^,...,v^} and 
set of lines E * {e^ e^jj}, the adjacency matrix A is defined by 


A * ' 1 n} , 

where a^j is defined as mentioned above. 

Clearly, there exists a one-to-one correspondence between the 
matrices A and F. Given A = w® set F^. * (a.^ .... ,a^.^}, 

F^j » {0,...,0}, if<j, and obtain F And conversely, given 

F = {F.j}, we have F.. as the i— row of A: 
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We note that the following normalized matrix has been 
Incorporated In the performance analysis: 


f|. • -7^ 


where 


and 


Ji- - 

’U 


0. If 1 j 


^1 “ Z ®1j* ^ • 

j-1 

This Induces the associated normalized adjacency matrix A. 


A - 


4^ *1j 

j-1 


But one observes that 


£ *1j ■ number of lines connected from node - deg(v^) . 
j*l 


We can thus write: 


“ 3igT^ 


11 


TT 


11 


deg(v. j 


We can distinguish between two main related performance analysis 
and design tasks. 

Task_A. Performance analysis and design of the antenna network, with 
respect to antenna radiation pattern characteristics (beam dispersion, 
etc.), when a proper topological configuration for the Interconnection 
of the power amplifiers Is assumed. 

Task B . The design of proper survivable families of topological 
structures for connecting the power amplifiers. 
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The study described here Is concerned with Task B. Incorporating 
the various Important topological measures presented In Section li.l 
we will develop survivable topologies satisfying a proper set of 
topological constraints. 

The network designer would then be able to choose the proper 
survivable cable connection configuration that satisfies the desired 
beam characteristics. Incorporating our Task B results, by Incorporating 
the latter Into his Task A performance analysis procedure. 

To Illustrate the significance of the synthesized topological 
structure. In yielding the proper measures of survivability, connectivity, 
number of lines, node degrees and network radius and diameter, we will 
present next certain optimal topological substructures. 

IL3 Survivable Network Structures with Minimum Number of Lines 
We have a set V of n nodes: 

V * Iv^ ,V£ , • . . , v^ } . 

We wish to connect them by a graph G that Is k-connected. There are 
many such graphs. As noted In Sec. II. 1 , each k-connected graph must 
yield a node degree now lower than k: 

deg(v^) ^ k , each 1 . 

Subsequently, the number of lines m(k,n) of a k-connected graph with n 
nodes must satisfy the relation: 

m(k,n) > 2 * 

This lower bound Is attained when all the nodes in the network have the 
same degree k. 


In fact, it can be proved that there exists a family of graphs, 
to be denoted as {H(n,k)), that attain this lower bound when nk is 
even. These graphs, attain the minimal possible number of lines of 
any n-node k-connected graph. This minimum number is given by 

m*(k,n) ■ [ J , 


where Lxj denotes the integer part of x. 

It can thus be seen that in a H(n,k) graph, either all nodes have 
degree k, or (n-1) nodes have degree k and one node has degree k+1. 

The H(n,k) graphs are constructed as follows. Let the vertices 
of the graph now be indexed as 0,1,..., n-1. Assume k > 2. 

For even k: 


1 1 if |i-j I ■ m (mod n), 1 m ^ k/2 
0, otherwise 


For odd k, even n : 


*ij 


1 , if |1-j I » m (mod n) , 1 ^ m ^ k-1/2 . 
or if |i-J[ » n/2 

0, otherwise 


For odd k, odd n: 


*ij 


1 » if |1-j I ■ m (mod n) , 1 £ m ^ k-1/2 
or if (i-j| » n-1/2 

0, otherwise 


We note that the adjacency matrix of an H(n,k) graph is 
constructed by following a very simple procedure. We have that A is 
symmetric, 

«ij * *ji 


and that 
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so that a-1 the diagonals of A are either all Ts or all 0‘s. 

Therefore, we need to specify only the first row of A In order that 
all of A (and the topological structure) be specified. However, the 
first row of A for H(n,k) Is specified simply by setting, synmetrlcally 
from both ends of the row, a total of k elements equal to 1. The 


elements inbetween are set equal to 0. 

Exanwle. n"7.k«4 



0 1 2 3 4 5 6 

0 

0 110 0 11 

1 

0 110 0 1 

2 

0 110 0 

3 

0 110 

A ■ 


4 

0 1 1 

5 

0 1 

6 

0 



Mote . Each vertex has degree k»4 

Graph has number of lirvs m«^nk*j* 7-4* 14 
Graph has diameter d > 2 


Graph is k»4-connected. There are 4 disjoint paths between 
any pair of nodes. 



- 142 - 





mtmi 



■1 



!n siMRury, network structures H(n.k) have n nodes, are k* 
connected, and as such have the minimal number of lines. If nk Is 
even, all the vertices In H(n,k) have the same degree k. If nk Is odd 


one vertex has degree l/tl, the rest have degree k. 
Example . n»7,k»3 



All nodes have degree k*3. except node 3 wMch has degree 
k+1-4. 

The graph diameter Is d*2. 



A 
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Graphs H(n,k) have a diameter which Is proportional to n/k. 

This can be large for large n. Thus, such structures will require 
long cables for our network, and will thus be Inappropriate. However, 
we will find that graphs H(n,k) are optimal for certain hierarchies 
within the global topological structure. 

To obtain a k-connected topology with minimal number of lines 
that also Incorporates radius constraints, as well as a set of central 
nodes, we next present the following optimal structures, recently 
developed and studied by us. 

II. 4 Reliable Centralized Topologies Under a Radius Constraint 

Assume now the nodes In the network to be decomposed into two 
groups. The first group is that of regular nodes, Vj^, the number of 
which Is n. We require each r3qular .“ode to have a degree not much 
higher than k, preferably not higher than k+2. 

The second group is that of central nodes, V^. Each central node 
is assumed to have a degree not higher than a prescribed number b, b > k. 
The number of central nodes Is set to be equal to s. 

We define the radius r = r(G) of the network G connecting these 
nodes, as the longest distance of any regular node from any central node. 
Thus, if the SPS antenna has radius r, then 

d^(v) ^ r, for each regular node v, 

where d^(v) denotes the distance from regular node v to any central node. 

We wish now to construct a graph which will be k-connected (reliable), 
have radius not larger than r, satisfy the node degree constraints, and 
have the minimal number of lines. 
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We have recently solved such problems, obtaining the following 
solution. For k _2, to obtain radius r and central nodes with 
maximal degree b, the number s of central nodes must be at least s^: 

* “ f BiT^TF)! 

where 

s(2,r) » r 

s(k,r) = ^ -^ ' 1 . k > 3 

It can then be shown that the topology with the minimum number of 
lines will then contain Sq central nodes, and have the following 
topological structure. 

We start from the Sg central nodes. Each of them Is then joined 
to the b distinct regular nodes, to form hierarchy 1 of nodes. Each 
of the latter Is then joined to (k-1) other distinct regular nodes to form 
hierarchy 2. We continue in this manner until hierarchy r is reached. 

The- resulting tree has a number of regular vertices not less than n. 

To obtain n regular vertices, we delete a sufficient number of groups of 
(k-1) end-nodes and add few more lines to obtain a graph where all 
regular nodes, except the end-nodes are of degree k. We have thus 
obtained a hierarchial tree T(n,r,k). The n^ end nodes of the tree are 
now joined as an H(ng,k-1) graph, which is (k-1 )-connected and has 
minimal number of lines. The central nodes are connected together as 
desired. The overall graph Is denoted as H(n,r,k) graph. We have 
shown It to virtually attain the minimal number of lines. Note that it 
is k-connected, has radius r and satisfies the nodal degree constraints. 
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The construction of an H(n,r,k) network is niustrated by 


the following example. 

Example . H(n,r,k) graph, with s-1 , single central node, r»5,k*3,b«4 



Note . #This network is k=3-connected. Between any pair of nodes 
there are at least 3 disjoint paths. 

•In this network, every regular node can be connected to 
the central node through a path which contains no more than 
5 lines. Network radius is r=5. 
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•Each regular node In this network has a degree equal to 
k"3, which is the minimal nodal degree for the network to 
be k»3-connected. 

•The central node has a degree b=4, which is its specified 
maximi degree. 

The number of lines in a H(n,r,k) network with s central 
nodes is noted to be given as 

m ^ j (k(n+l) + j s(b-k) + , 

where m^ represents the number of lines used to connect the s centers. 

These H(n,r,k) topological structures and various proper 

extensions of them to be developed in our topological design studies, 

will serve as fundamental hierarchial topologies in the construccion of 

the desired families of survivable topologies. 

1 1.5 Survivable Topological Structures for the Antenna Phase-Array 
Network 

In arriving at the optimal topology for connecting the antenna 
center to the distributed power amplifiers, we are constrained by the 
existing locations ofthe amplifiers (regular nodes) and the center. 

Thus, we are bound by the geographical locations of the network nodes. 

Subsequently, we are able to classify the network nodes into r 
levels, so that every node at the i-th level is at a distance of i 
distance units (or approximately so) from the center. Connections will 
generally be made in such a manner that usually a single cable piece 
(line in the underlying graph) will cover a single distance unit. On 
certain occasions, a cable piece will extend over a few distance units. 
Subsequently, a graph structure with minimum number of lines, satisfying 
the above distance structure constraints, will generally attain also the 



minimal overall cable length. 

To Illustrate the geographical constraints, we draw again the 
H(n,r,k) network In the example In Section IV, Incorporating now the 



Recall that the Illustrated structure is 3-connected; each node 
has degree 3 except the center, the degree of which is 4; it has radius 
5; and it contains the minimum number of lines among all such graphs. 

Of course, in the actual tree network structure the number of 
nodes in each level is already prescribed. This will be taken into 
consideration in developing the optimal structure. 

The H(n,r,k) structure is k-connected, so that upon the failure 
of (k-1) or less nodes every remaining pair of nodes is able to 
communicate; and in particular, every remaining node stays connected 
to the center. This network has radius r, so that every node has a 
path not longer than r lines which connects it to the center, under 
no failures. However, under any failure, even though a remaining 
node stays connected to the center, its distance from the center will 
usually increase. This may induce an undesirable extra phase offset 
at such a node. 

For example, node V 2 is at distance 2 from the center (through 
vertex v^). If, however, vertex v^ fails, node Vg will be at a 
distance 9 from the centerl Thus, it will have to experience a 
considerable possible phase offset due to the additional 7 lines in 
its shortest path from the center. Similarly, the distance of v^ from the 
center increases from 3 to 8, upon the failure of 

To prevent excessive additional phase offsets due to node failures, 
we need to develop topologies which are radius stable. Such structures 
will guarantee the distance for many node to the center to remain the 
same, or increase only ina restricted manner, upon certain node failures. 
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We have recently studied the topological design of radius stable 
reliable topologies. Results of these studies will be coupled with 
the considerations and results mentioned above to develop optimal 
tree network structures. 

For example, It Is possible to design topologies where certain 
node level assume proper H(n,k) structure, while the levels are 
connected approximately according to an H(n,r,k) topology. Such 
structures can be shown to have certain radius-stability properties. 

For example, the following figure Illustrates a 3-connected 
network structure which has certain radius stability properties. The 
network does not necessarily contain the minimum possible number of 
lines. Note that upon the failure of a node or two, the distance of a 
node from the center does not usually increase considerably. 

II. 6 Recommendations 

It Is our objective. to Incorporate the powerful methods, techniques 
and results mentioned here to derive optimal survivable topological tree 
structures for connecting the SPS phase control centers to the central 
(master) node and for connecting the power amplifiers. 

In deriving these topologies we will impose the following topological 
constraints, which we have shown to be closely related to the desired 
phase control system performance objectives. 

• Tree (deterministic and probabilistic) connectivity and 

Invulnerability measures. 

•Nodal degrees. 

•Number of network lines and overall cable length measures. 

•Network radius constraints and radius-stability measures. 

We recommend subsequently developing families of optimal topologies 
satisfying the above prescribed measures. One would then be able to 
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incorporate these proper survivable structures in the perfoiTnance 
analysis of the reference phase distribution tree, to conclude the 
associated beam dispersion characteristics. The topology yielding 
the desired dispersion properties could then be chosen. As such, 
our study will guarantee the latter to be highly reliable, highly 
insensitive to node failures and requiring the least amount of overall 
cable length and the least number of nodal connections. 
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APPENDIX III 

NSTRS: A NEW METHOD OF TIME (PHASE) 
TRANSFER FOR THE SPS PHASE CONTROL SYSTEM 



III.l Introduction 

Most tine transfer systems In use today utilize MS (master to slave) mode 
of operation mainly because the alternate ways of achieving the 
time transfer such as Mutual Synchronization methods are not fully 
understood. Recently, there have been some attempts made to under- 
stand the mechanization of time transfer In a mutually synchronized 
system Refs. [1-3]. The steady state frequency and phase errors 
of the nutually synchronized systems are, as pointed out In Ref. [3]. 
very sensitive to the path delays in the system, hence the necessity 
of delay compensation arises. Recently, the delay condensation Is 
attempted In Refs. [4 - 5]. 

The Returnable Timing System (RTS) described In Ref. [5] Is of 
particular Interest because the delay compensation feature of this 
system can be combined with the MS system to create Master Slave 
Returnable Timing System (MSRTS), described In this paper to combine 
the advantages of MS as well as RTS configurations. MSRTS uses the 
basic MS hierarchlal structure (tree structure) to transfer time 
from one station to any other station perfectly (delay Independent) 

In a noise free environment. This method Is particularly Interesting 
because by addition of some equipment, the already existing MS systems 
can be made delay condensating and the theory developed here can be 



ustd to prodict ptrformanct of th# systom. 

Tht basic Idea In MSRTS Is to send the signal received at a 
station via a delay back to the sending station via the same delay 
(path). The delay accumulated by this returned signal Is used to 
advance the phase of the master (sending) station thereby canceling 
the effect of the delay Introduced by the path. To understand this 
Idea clearly, we shall start with only two stations one of them Is the 
master station and the other one Is a slave station. 

III. 2 Mathematical Model 

Fig. III. 1' depicts the system of two stations In a chain of stations 
completely. Each station has a signal processing equipment with a 
phase locked loop In It. Since our paper Is primarily concerned with 
the synchronization of the stations, we will represent the stations by 
their respective phase locked loops only. 

As described In the figure, we desire to transmit the phase of the 
master signal to the slave station. The master signal and the returned 
signal from the slave station are used to drive the PLL and the output 
of the VCO Is sent to the slave station via a directional coupler. 

This directional coupler has the ability to sense the direction of 
energy flow and hence It can route the signals coming to It at the 
same frequency In the pioper directiws. At the far end of the cable 
there Is a directional coupler with termination. The termination In 
the coupler simply reflects some of the Incident energy back to the 
master station for the purpose of delay condensation. It should be 
noted that the mechanization can be brought about by using two 
diplexers (or four bandpass filters) Instead of the two directional 














Master Signal 


s(t) » A sin(o)t+ 0 |^|) 


(III.l)' 


where is the phase to be transferred to the slave state 

r(t) » K cos (ut + 0 (t)) 

where e(t) is the estimate of the master PLL of the phase 6 u 
s^(t) = A sin 2((Dt+0|^|) 

S 2 (t) = AK sin{o)t + 20^^ - 0 (t)) 

SjCt) » K cos(u(t-T^2’^21^ 

where t-j is the delay suffered by the signal from the j— station to the 
1 — station 

c(t) = AK sin(20n^-0(t-T^2’'^21^"®^^^ * "^■^12‘^‘"21^^ 

and with no frequency offsets and no biases, in the steady state we will 
have 

0(t-Ti2-T2i) - 035 

° * 2®M " ^^®ss ■ “s^l2 ■ “s^21^ 

* 20^ - 2635 + U5T12 + (O3T21) ; mod 2 tt 

where 033 is the steady state value of e(t) and is the steady state 
frequency. 

®SS ' °H *"S — V-- ” U”- 3) 

Thus the phase ambiguity goes away when multiplication is an even 
integer. 




steady 

state 


r(t)l 


» • K cos(a,5t^^^) 


K C0S(a)^t+6|^^ +u^ 


Phase of at point B (at the slave station) 


, + ru*^2}\ 

“s V 2 7 ’V21 

“s y 2 “ ■^2lj 

, /^12’‘'2l\ 

M ^“s\^— rH 


= 8 


“ «M " “s ‘t> 


(III. 4) 


where 


At ^ t ^2 ■ ■*^21 


(III. 5) 


Now At=0 because the same signal travels back and forth through the 
same line, therefore ti 2 ®t 2 i* (This assumes a nondispersive cable). 
When at=0, in the steady state, the phase of the master signal at 
point A equals the phase of the received signal at the slave station 
at point B. 

III. 3 Chain of Stations 

As shown above, the phase of the master signal is regenerated 
at the slave station (point B in FigJILl. This signal will be fed to 
the slave PLL and this ’’LL will act as the master for the next station 
in line. FigJII.2 shows the situation in a block diagram form. Phase 
of the signal at the input of any slave (in the steady state) is 
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obtained from the phase of Its immediate master by properly Interpreting 
the quantities In eq.III,4. For example the phase of the signal Input 
to the (n+1)— slave In the link Is given by 

Input Phase to . . ^^^12 ^^‘^23 ^ ^^’'n-l.n ^ ^^n,n+l 

the n+lst Station ®M ~2" 2 '* 2 2 

1=1 

where Is the phase of the master oscillator to be transmitted to the 
slave stations. Ideally, with nondlspersive cables, we will have 

i +1 *■'1 i+i ■‘*1+1 1 * ® VI = 1, ..., n»l, hence (III. 6) gives us 
the phase Input to the (n+1)— station Is the same as the master station, 
III, 4 Imperfect Compensation 

The condition * 0 for the Ideal case assumes not only the 

cable be nondlspersive but also that all the components In the system, 
especially the directional couplers are perfect. In actual practice, 
this Is rarely the case hence It becomes reasonable to assume that each 

•At. is a random variable with zero mean (unless the system Imple- 

’ 2 
mentation dictates otherwise) and a variance, say It is also 

reasonable to assume that there Is no correlation between any pair of 

random variables. In fact, we will assume that they are statistically 

pairwise Independent random variables. Thus 

e[ati,i^i] = 0 

^^^■^1,1+1^^ ' ®i,i+l 

with 

E{[AT^^,.^^][ATj * 0 Ifj (IIIJ) 


2 

®i ,1+1 


and using (III. 7) get 

E[yJ] » Var[y] = \ 

1*1 



^ o' 

y 


(III. 9) 


n 


Thus the disturbance to the phase Input to the (n+1)— station has 


a zero mean and a variance of o' . A simple application of the central 


• I 

limit theorem gives that for a large enough ''n'' we have 


pw (y) 


7^ 


-y'/'oJn 


(III. 10) 


and 


Phase Input to 


the (n+l)s^ Station “ *n+l * ®M ^ ^n 


(III.ll) 


then 

ECVl’ ■ ■ 'n 

and 

Var[« _,] “ o' (III. 13) 

n 

111,5 Effect of Nonzero Mean of At ^. on the System 

So far we have tacitly assumed that = 0 Vi ' l,...,n. 

If there Is some reason to assume that E[At^. f 0 then 

the mean of y increases In the following way; 
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E[y„] 


i'n ■ 

i-1 

?Z “ ?E *^1.1+1 ° 

1-1 i-1 


- ifc 


1,1+1 


c Vi. 


- I E <'1,1+1 

1-1 


(III. 14) 


In such a case another mechanization of the system described In F1gJII.3 
may prove more useful . 

The key points to note In this mechanization Is that the directional 
couplers In Fig. IIL3 are replaced by diplexers and frequency dividers are 
Included. Also the direction of flow of "u" signals and ''2(u" signals 
Is alternated. Going through a similar analysis as before we can arrive 
at an equation similar to (III. 11) giving us the phase of the input signal 

St 

to (n+1) — station In the steady state to be 







c Vi ■ then 


'l.i+l 


ECy„3 ■ H)'*'. 

i-1 


0 if n is even 


if n is odd 


% 


i«l ” 

Con^aring (III. 14} with (III. 17), Me see that the mean of the phase disturbance 
has decreased by at least a factor of n while the variance remains 
the same. 

III. 6 Probability Distribution of the Phase 


input phase to , 
(n+1 )U Station 

■ Vl • 


where 



^n 

with 

1 A 

“ rE^^i.i+l 

i-1 


i-1 

i.Hl^ - ? t 

i-1 

S.i+1 ^ <=n 

and 



V»r[y„] 

1 2 
T L ®i,i+l 

i-1 

^n 

EtV)] 

• 9„ + ECy„] - 

^n 


n 


Assuming a large enough 

"n" using the central limit theorem, we 






(III.18) 
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Thus It becomes possible to find the probability that the phase Input 

kSt 


to the (n+1)— station lies between +a of the master signal phase, e^. 
This Is given by 


1 9n+l - 




'/j 


'M 


{ 


M 


,+A 




7 “ " V V 




^n+1 


erf 




» erf(a-c /o„ ) + erf(A+ c„/o„ ) 


where 


n y, 


erf(x) 


- yLf 

72P 

n 


X 2 " 

e-y 


(III. 19) 


Thus we see from (HI. 19) how the mean value c„ enters into the equation 
to affect the probability of being In (-a,a) interval. If c^»0 then 
(III. 19) reduces to 2 erfCi/oy^). In such a case with P[6n^-A ^ 

Vl ‘ * \ ’ i ‘’i.i'+l * ® 

we can get from (III. 19) that 


a - 2 erf 


-jO. 


n « 4 X (A/a)^ X 




(III.20) 


The above equation tells us how many stations ("n") can be allowed in 
a chain of stations for fixed values of "a" and "o" and a selected 
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value (variable) of a. The ratio {Uay can be intrepreted as 


(J)' ■ ‘ Spr... <”>•*'> 


If we let A>ir and change a to (1-a) (20) becomes 


n 


4( jr 




(IIL22) 


“n" gives us how many stations can be allowed in a chain for a pre- 
determined value of "a" and the system parameter value o before slip 
of one cycle is expected to occur at the input to (n+1) station. Thus 
a large number of stations can be allowed if o is sufficiently small 
before the cycle slipping occurs, alternatively the probability of 
slipping a cycle tends to zero. 

III. 7 Phase Noise Propagation in a Chain of Stations 

The above analysis does not consider any form of noise either 
generated by the system components like the PLL or the channel noise. 

The analysis presented below considers the propagation of noise generated 
by the Pll through a chain of stations. Fig.III .4 describes the connection 
and signals at the k~ st tion. We will adopt a new notation 0) means 
phase of the signal (•)• At the k— station 


rjj ■ ut + &|j(t) 

where 0j^(t) ^ 9 |^(t) + ♦,^(t) in which ^j^(t) is the phase noise generated 
by the k— osctllator and 0|^(t) is the phase estimate of the k— 
oscillator. 
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Figure III. 4. Connections at the k— Station in a Chain of Stations 
















\{t) - 

H (t) - «t - 2 wx,j c IJ.,) - e|((t) 

1 

1 (t) • «(tt|,*,,|, - Xk.i,*,) ♦9kf‘-’|i*1.k-'k.lc*1> 




Taking Fourier transform on both sides, 

• ^“(Tk+i.it+'k.k+r^’k.k-l^ * ^ ’®k^^ 

./i-^^.l.k^k.k-l^ -e^L.) 




*k,k-1 




Vl,k^%,k+l^j 


kF. (tt) 

0j^(«) ■ ju" Cc|^(m)] »k(“) 


kF.(u) , . , , k /'^“■^k.k-l 

®k^“^ ’ ”17“ t^’’“^Vl.k*^k.k+r^’k.k-r ^©k.iUJe 

-e^U) (I *e‘^“‘'’‘*’'’'‘ ’"■'’)] *»k'“> 

Assuning that proper filtering removes thes’s and solving for 0 ^(„) 
from ( in.23) yields 


kF^(o,)e 


'k.k -1 




®k-l^“^ 




“ ' IuTtT k+l^XL+1 kS 

kF ''’ ’ 3 (III. 24) 



Ju ♦ kF|^(u)[l+e k+i»k k,k^l ^ 


H.(«) 


.Ju(- ^ p- (III. 25 

J« ♦ kFj^(«)[l+e ^ k.k+l^j 


Then (III. 24) becomes 


ek(«) • R,j(«)0^_,(w) ♦ H,^(u,)<»^(w) (III.26] 


This Is a recursive relation for e,^(«) and can be solved for the n^ 
statloti phase output. This turns out to be 

■ {S V“) * i {"((“) »j(»)|*,(.)(im7) 


where 


Thus the system can be represented as In Fig. m.5 in a series form. 
Or this can be modelled as a parallel fed system shown In Fig. III. 6 
From (II I. 27) we have 

®n^“^ ■|k|' ^ ®n^“^ • E ^T^^RjU)|<-^(a,)(III. 28 ) 

Now assuming that E[a^(a.)] - 0 Vi and E[»^(a,)^j(u,)] - 0 ifj we have 


£[♦-(«)] ■ o»^e_(«) » TT r.Il.1 qJu>\ 


1111 
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Figure III. 6. Parallel Equivalent of the Chain of Station* Connected in MSRTS 










i.e. OP the average, the output phase of the n— station oscillator 

equals the master phase. We would be interested to find o . 

♦n 


S (w ) = 
♦n 




where S (u) is the single sided spectral density of the ij).. 




■ 27 


' 27 f (u.)da, 


f Z|H,(a.)|2 ! fr |R,(<u)|^is. U) 


I !Hi{fi>)| S (wjdu 
Jn ^ 


where R|^(u) and H|^(u) are defined in 011.25). How if we assume that 
T^j= T Vi , j we have (u) = H(u) Vi and R^(u) = R((d) Vi. In such a 


r |H(o))|^[|R 1^^'’'^^S, + |R|^^""^^S. +... + |Rl^S, 


lHU)rs^{u.)d 


How we let S, (a») = S,(oa) Vi hence 


|H(<o)r S.,.(u;)d 


Z 
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Now assuming T^j»0 VI J and F|^(u) » F(u) VR and S^(u) = 1/u we 
get 




“ JOi ■+“If{o)T ’ 


Assuming 


FU) = T-+ 


1 + juT. 


1 + j(i)T 


1 


R{-) 


-■? where x = o)/u_ 

(1-x^) + j2;x " 


with = /R/ T^'* and c = gives 

2 

H(ai) = 1 - R(u>) » ^ 

(1-x^) + j2i;x 

Substituting H(a)) and R(u) in (III. 33) one can compute the normalized 
mean squared random jitter after n stations (normalized by the 
mean squared random jitter produced by one station) is computed by 
the use of an IBM Computer. The re'^^ults are presented in Fiq. III. 7. As 
can be seen, as c, the damping factor of the circuit increases, the 
noise variance at the end of n station chain decraases but this may 
create problems in acquisition mode of the stations. 

II1.8 MSRTS Tree Network 

When there are a number of stations to be connected (phase locked) 
together to transmit phase information from one station to all the 
remaining stations, one possible way of connection is to connect all 
of them in a chain and let the master feed the chain. There are 
obvious disadvantages of a series connection with the above scheme. 
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hence the well known master-slave hierarchial tree structure will 
be tried next. FigJII.8 shows such a tree structure using MSRTS 
concept for the time transfer. 

The master (station 1) in the figure feeds the phase to four 
physically close equidistant (from the feeding point) phase locked 
loops via a four way power splitter. These are level 1 slaves. Each 
of the level 1 slaves is connected to a four way power splitter as 
level 2 via the MSRTS. Again a group of 2— level four PLLs are 
placed equidistant from the power splitter and each of these is 
connected to the third level power splitter via MSRTS. This is 
continued until all the oscillators (stations) are connected. 

For each oscillator in the network there exists only one chain 
connecting the master to the oscillator and our analysis in the 
previous section holds. 

Ill .9 Conclusion 

MSRTS provides a new way of delay compensation which can use 
already existing Master Slave systems. Proper design of filters at 
each station can reduce the varianceof noise at the end of the chain 
of stations. MSRTS can also be useful in antenna theory, especially 
those reflector antennas where distribution of the master phase 
accurately is vital to the performance of the antenna. Also this 
method can be of use for sensing phase at a remote location for phase 
measurement purposes. 
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APPENDIX IV 


SOLARS IM SOFTWARE DEVELOPMENT 

This section describes the SOLARSIM programs and SOLARSIM 
capabilities to quantify the spacetenna performance parameter values. 

The overall setup of SOLARSIM computer simulation Is shown In Figure 
IV. 1 where the user Inputs direct the central processor unit to pick 
the correct SOLARSIM subroutine from the subroutine package stored In 
the memory of the computer and after computation, produce the desired 
output either In the form of plots or printout or both. Figure IV. 2 
shows the capabllltels fo the SOLARSIM subroutine package. As shown 
In the figure, SOLARSIM can compute the following quantities; 

1. Antenna element covariance matrix. 

2. Individual realization of random power patternT 

3. Mean far field power pattern. 

4. Main beam gain losses. 

5. RMS pointing error. 

6. Tilt/mechanical error effects on gain with effect of 
conjugation at other than radiating element level. 

7. Main beam power transfer efficiency. 

There are two different types of inputs necessary for ^he operation of 
SOLARSIM subroutines, users inputs and the computer generated inputs. 
Both these Inputs are shown in Figure IV. 2. The computer generated 
inputs are automatic in the sense that once the user's input specifies 
the SOLARSIM subroutine, the subroutine calls for the necessary computer 
generated Input for its execution. 
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USER'S INPUTS 

PILOT SIGNAL ORIENTATIONj 

PHASE CONTROL TREE 
STRUCTURE 

•PHASE JITTER STATISTICS 

•AMPLITUDE JITTER 
STATISTICS 

•CURRENT TAPER 

•POWER TRANSPONDER 
CHARACTERISTICS 

•MECHANICAL STATISTICS 

•MEDIA 


r.nMP[ITER riFNFRATED INPUL 

•LIBRARY SUBROUTINES 
•PLOTTER INTERFACE SOFTWARE 
•PSEUDO RANDOM NUMBERS 


SOI ARSIH ROUTINES 



•ANTENNA ELEMENT COVARIANCE MATRICES 

(aecm) 

•RANDOM POWER PATTERN SIMULATION 

(rpps) 

•EXPECTED POKER PATTERN (EPP> 

•MAIN BEAM GAIN LOSSES (MBGL) 

•RMS POINTING ERROR (RPE) 
•tilt/mechanical ERROR (TME) 

•MAIN BEAM POWER TRANSFER EFFICIENCY 

(mpte) 


Figure IV. 2. SOLARSIM Capabilities. 
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The SOLARSIM programs are configured such that all the user has 
to know are the Inputs to the programs. In certain programs this step 
is made easier by having a question and answer session in the beginning 
of the programs. 

The documentatiwi below contains information necessary for the user 
to use the programs. The topics covered can be grouped into three 
headings. The first part gives an overall input output diagram for 
the program along with a short description of what the program does. 

The second part shows the flow chart of the program and the third part 
shows a typical output of the program. 
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SOLARS IH Subroutine PVRPAT 


Purpose ; PVRPAT slimilates the SPSS antenna system and computes the 
actual (not the averaged) power patterns for the antenna, when 
the phase distribution system noises are present and not present. 
This program gives a probable pattern for the SPSS antenna. 


Inputs ; 

L - 


LL • 

OB - 
DX > 

DY - 
STEP - 

NN • 

THETAS - 
THETAO - 
PH 10 - 
PHI - 

SGMA - 

SEED - 


Number of levels of the phase distribution tree 
structure. 

Nim^er of current levels. 

Current taper In dB. 

Incrmnental distance In x direction. 

Incremental distance in y direction. 

Increment in 6 direction for the computation of 
power pattern. 

Number of d values to be used to compute the 
power pattern. 

Starting value of e direction 
Pilot signal 6 direction. 

Pilot signal ^ direction 

Constant value of ^ direction to be used In the 
power pattern. 

Array for the standard deviations In each branch 
of the phasing tree (In degrees). 

Array of Initial values for generating the Giusslan 
random variable representing the phasing system noise. 




om 


Outputs ; 

T - Array of e valuas (minutes). 

XC- Array of u ■ sin e cos ^ values 

YC> Array of computed power pattern. 

A simple output of the PVRPAT Is shown In Figure IV. 4 where 

L ■ 6, LL ■ 4 and dB ■ 10. There are two columns under "ratio". 

The first one Is the antenna power pattern with no phase distribution 
system error and the second one Is with the phase distribution system 
errors specified In the Input. Shifting of the maximum of the main 
lobe Is evident In the second colwwi of “.‘atlo". This shifting of 
the main lobe Is plotted In Figure IV. 5 for the case of L ■ 5, 10 dB 

current taper and the phase error standard deviation of 10* at each 

level of the phase distribution tree In this particular case ihe main 
lobe has shifted from Its mean position by 0.3'. 
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INPUTS 


k > 0 


INITIALIZATION & SCALING 


SET UP SOUNOARIES 
FOR THE CURRENT TAPER 


RANDOMLY GENERATE A PHASE JIHER 
AT THE END VERTEX OF THE TREE 


DECODE THE CONNECT POSITION 
OF THAT VERTEX IN THE ANTENNA 


DETERMINE THE CURRENT 
MAGNITi.’^t AT THAT POINT 


COMPUTE THE POWER 
PATTERN 


k < A"- ? 


repeat for DIFFERENT DIRECTION 


Figure IV. 3. Flow Chart for “PVRPAT" Program. 

- 187 - 


JZinCZ> 






. »>}* i 


t *# i *( I *# l * r !* fl * 

i«it.iii/iiu«iiUi ,miii 
) j > r.xii 

LlhKU»3 PiEU»MA 

> 


eXCT 

iLM'24f<2 «<L7l-3 ••7/17/7H 11:36:52 
i:TAPY»»>riit33, PrtOC SIZE (1/U) *5690/11009 
3Y.:5*i LI0V. LEVEL 72-C 
GKO .-SAP 


HiAUb bATUrl 


TflETAI' PlliG RO SCMA(i) DO DX 

.i?i:0 .000 320.000 l.OOe 10. 00b 1.000 


311ETA U 


RATIO 


.CEOdCrOO 

l.ocooeool 

2.00000003 

3.00000004 

4.00000005 


.Kooooeofj .noroncB 

.01492625 IO37091+COG 

.G2985251 -.4157336+000 

.04477676 -.9373706+000 

.05970500 -.1670514+001 


921421. J001C014 3C16S9. 95946666 

Figure IV. 4. A Sample Outpuf of PVRPAT. 



oy LL L 
1.000 4 6 


-.5747175-001 

.0000000 

.1494701+000 

-.5096344+000 

.1086073+001 
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SOLARSIH Subroutine EAPP 
Purpose ; EAPP computes the exact expected (averaged) power pattern 
for the SPS antenna for different configurations of the phase 
distribution system. It also Allows the current taper to be a 
variable of the program. 

Inputs ; Figure IV. 6 describes the inputs to the EAPP which are 

i 

basically the same as the list of inputs to the PVRPAT except 
for the last entry which is changed to SGMD. Fig. IV. 7 gives 
the flowchart of EAPP. 

Outputs : The list of outputs is similar to that of the PVRPAT output 
list. 

A sample program is run for the case L = 5, LL = 4, dB = 10 dB. The 
output e is in minutes, u » sin e cos <|i and the last column gives the 
normalized gain of the antenna. Figure IV. 8 gives the particulars 
of the sample run. 










USER'S INPUT 


►ANTENNA DIMENSIONS 


•ANTENNA ELEMENT POSITIONS 


•PILOT SIGNAL ORIENTATIONS 

•CURRENT TAPER AND LEVELS 

•PHASE DISTRIBUTION NETWORK 
STRUCTURE 

•PHASE JITTER STATISTICS 
(cam be DIFFERENT FOR EACH 

branch) 



Figure IV. 6. Inputs to EAPP Subroutine. 
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#4#10#O» fSflti* #0* fO« #!• #320* 

>SwS*.0*»0-#S. 


OWO'NAL PAGE » 

OF POOR QUALITY 


ltCMA( 1 ) • 
SXGMA( 2 ) - 
SICHA( 3 ) - 
8ZCKA( 4 ) - 
SICMA( S ) > 


s.eeeeeoa* 

s.eeoeeooe 

s.eaaaeoaa 

• 

s.aaaoooea 

s.aaaaaaaa 


THETAO PHia RO DB DX OY LL L 

.eaa .aae 328. oea le.aae i.aeo i.aaa 4 s 


THETA U 


(PERTURB) RATIO 


.aBaaeoeo 

le.oaaooooe 

ae.aaeooooo 

29.99999976 

39.999999S2 

49.99999905 

59.99999905 

69.99999809 

79.99999809 
89.99999689 


.aaeeeaee 

.14926883 

.29852040 

.44777745 

.59783071 

.7462789B 

.89552078 

1.04475510 

1.19398056 

1.34319593 


.0000000 

-.4152529401 

-.2011022402 

-.1548799402 

-.2231569402 

-.2111276402 

-.2155554402 

-.2391221402 

-.2326496402 

-.2579810402 


364244.17187500 


Figure IV. 8. A Sample Output of EAPP 
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SOLARSIM Subroutine APPRXPAT 

Purpose ; The previous subroutine computes the actual averaged power 
pattern but It Is slow and uses a lot of CPU time. APPRXPAT 
computes the averaged power pattern by Introducing an approxi- 
mation to the phase error covariance matrix which makes the 
computations considerably faster and the results are very close 
to EAPP results. APPRXPAT, like the EAPP, handles various configura- 
tions of the phase distribution system and also different current 
tapers can be specified. 

Inputs : 

XRM - Alphanumeric array for generating the format of the 
output printer plot. 

L - Number of levels of the phase distribution tree 

structure. 

LL - Number of current levels. 

NTER - Number of 0 values to be used to compute the power 
pattern . 

NB - Array containing number of branches per level of the 
phase distribution tree structure. 

SGM - Array of standard deviations In each branch of the 
phasing tree (degrees). 

START - Starting value of the 0 direction. 

PHIO - Pilot signal direction. 

THEATAO- Pilot signal 0 direction. 

PHI - Constant value of direction to be used In computations. 
Current Taper In dB. 


DB - 


Figure IV. 9 gives the flowchart of the subrouting APPRXPAT. 
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READ IN DATA 


ORIGINAL PAGE IS 
OP POOR QUALITY 



Figure IV.S. The Flowchart for the APPRXPAT. 








ns^om 


Outputs ; 

T - Array of the theta values (minutes). 

EFF - Array of the approximate power pattern. 

Library Subroutines: PLOTl 

Figure IV. 9 shows the sample output run for the case of a four 
level phase distribution tree with 10 dB current taper Induced by 
10 current levels and no phase errors (I.e., L>4, LL ■ 10, dB ■ 10 and 
sigmas ■ 0). The branching of the phase distribution tree Is 
16 X 16 X 16 X 25. The last two columns represent e and the associated 
approximate normalized power pattern. 

Figure IV. 10 depicts the plotted output of the APPRXPAT for the 
four level phase distribution tree structure described above with the 
standard deviations of phase error at each level equal to 7.5®. The 
current taper Is 0 dB for the solid line curve and 10 dB for the dashed 
line curve. As can be seen from the figure, there Is a widening of the 
half main beam width from 0.58 minute to 0.67 minute and a drop in the 
first sidelobe level from -13.4 dB to -17.5 dB. 

Figure IV. 11 shows the power patterns for a nine level phase 
distribution structure with four branches per node and 0 dB current 
taper. The standard deviation of the phase error added at each of the 
nine levels Is a variable. The values are given in the figure. As Is 
expected from the theory, we see frcmi that the main lobe becomes wider, 
loss of gain Is experienced and the nulls begin to fill up as the phase 
error variance is increased. Number 5 curve with 90® standard deviation 
of phase error at each level shows that the antenna Is no longer 
distinguishable from an omni antenna. 
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Figure IV. 10. Power Pattern by APPRXPAT 









SOLARSIH Subroutint TkUEGGO 


Purpoie ; TRUE6G0 computes the actual main beam gain of the SPSS 
antenna with and without the noise Introduced by the phase 
distribution system. 


Inputs: 


L - 

Number of levels of the phase distribution tree 
structure. 

LL - 

Number of current levels. 

OX - 

Interelemental spacing In x direction. 

OY - 

Interelemental spacing In y direction. 

S6MA1 - 

Array for standard devistlons in each level of the 
phase dsitribution tree structure. 


Outputs ; The output consists of two gains, one when the phasing 
system Imperfections are present and the other when there are 
no Imperfections. Ratio of these two generates the normalized 
gain sample run shows the output for L ■ 4, LL ■ fi and o for all 
levels • 5* 

Figure IV. 12 Is the flowchart of the TRUEGGO program. Figure IV. 13 is 
a plot of normalized gain in dB and the total phase error at the 
tip of the tree (equal amount of phase error added per level ) as a 
result of the sample run of the APPRXPAT. The curves are drawn for the 
nine as well as four level phase distribution structure with current 
taper of 0 and 10 dB. As can be seen from the figure to achieve a 
gain efficiency of not less than 96S the total phase error at the tip 
of the tree cannot be allowed to becor« more than 12”. 

APPRXPAT subroutine was used to determine the effect of the phase 
error variance added by the last level of the phase distribution tree 


structure. A sample run of APPRXPAT for the four and nine level phase 
distribution tree structure was made with the last level sigma varying 
^ (the sum total of remaining level sigmas was held constant to 15®) for 
current taper equal to 0 and 10 dB. The resultant plot Is shown In 
Figure 'IV. 14. 

Effect of phase error added by different levels In the phase 
distribution tree Is shown In Figure IV. 15. Two different structures 
are considered. One with eight levels of phase distribution system 
and the other with nine levels. The normalized gain Is plotted when 
the standard deviation of the phase error added by the system at the 1— 
level is 5° while the remaining levels add no phase error 1 = 1,...,9. 

It can be deduced from this figure that as far as the gain reduction is 
concerned, the last level is the most important level, i.e., t.he phase 
variance of the PAs will have a worse effect than the oscillator noises 
added at the first few levels. Fig. IV. 16 is an extension of the above 
idea applied to a nine level phase distribution tree. Curve 1 is a 
plot of normalized gain with respect to the first level sigma varying 
(all other levels add no phase variance). Curve 2 is a plot of 
normalized gain against the RMS phase error build up with each level 
sigma varying. Curve 3 Is the plot of normalized gain against phase 
error added at the ninth level (all other levels add no variance). 
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SOLARS IM Subroutine APPDEL 


Purpose : Given the phase distribution system, the imperfections 
(noises) associated with it, APPDEL computes the value of 
standard deviation of the pointing error. This program also 
takes into account the specified current taper. Figure IV. 17 
shows the inputs necessary and the outputs generated by APPDEL. 

Inputs : 


L - 

Number of levele of the phase distribution tree 
structure. 

LL - 

Number of current taper. 

DB • 

Current taper in dB. 

LB - 

Array for number of branches per level of the phase 
distribution tree. 

ICNTL - 

Control parameter. 

1. For printout with supplied sigmas for the distributi 
tree. 

2. For plot of (S.D. of pointing error) vs RMS . 

3. Stops the execution. 

SGMS - 

Array for the standard deviation in each branch of 
the phasing tree. 

NTER - 

Number of points desired on the plot. 

ICNT - 

Control parameter. 

1. For plot of vs first level . 

Outputs: 

2. For plot of vs RMS a with all levels varying. 

SIGLOW - 

Array for the lower bound on the a. . 

SIGDLU - 

Array for approximate 


ISGDUP - Array for the upper bound on 


Library Subroutines: Plot 1. 

The program has question and answer session, the output Is self 
explalnatory. Fig. IV. 18 Is the flowchart of subroutine APPDEL. A 
sample output Is shown In Figure IV. 19. Figure IV. 20 Is a plot of 
the RMS pointing error against the total phase error obtained from a 
typical run of APPDEL for four and nine level distribution trees for 0 
and 10 dB current tapers. The taper broadens the beam but It becomes 
more susceptable to pointing error also this can be ascertained from 
the two curves In Figure IV. 20. As the figure shows for 10 dB 
current taper for either four or nine level phase distribution 
structure for an allowable probability of 0.95 for the absolute pointing 
not to exceed 4% of the 3 dB beamwidth the maximum phase error build 
up can equal 10®. 

Figure IV;20 shows the plot from another sample run of APPDEL. 

The phase distribution tree structure has eight and nine levels only 
one level had a phase error input of 5® all the remaining levels added 
no phase errors. The importance of the phase errors in the first few 
levels is evident from the figure. 
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Figure IV. 18. Flowchar^^of Subroutine APPDEL. 



















lO.tSftO 

>4»4*t4tl4*29 


10.00000000 


,99427&59 .9/177283 .92828106 .86666746 .79082099 .70929383 .61477319 .53374151 .43609009 


«10799472«07 


. 10606070109 


•5S74090410# 


TYf^E 1 FOR FOfNTOUT WITH SUm.lCD SlCnASt TY9C 2 FOR FtOT OF ICtU VS THC ktlS SIOHO* TY9C 3 TO ST09 

>1 

TVf'E IN TtC SICnAS 


.boioWfV 

.02925606 

.07809840 

1 FOR FRlNrOUT WITH SUPf-LlCP SIGHASt TTFC 3 FOR PLOT OF »CLU VS TNC MtS SlGMR* TYPC 3 TO STOP 

IN THE HUnl>eR OF l-QINTS tCSIRCD ON TiC PLOT 

IN THC INCRCnCNT FOR X AXIS 

I TOR F*LOT OF ICiU VS 1HC FIRST LFVCL SICNA OR 2 FOR I«« U VS RNS SIOnR VITH ALL LEVELS VARTIND 


1 .OOOOO O 4 

2.00000 I l« 
l.'OooOO T 

OOOOO 1 

5 . 00000 I 
6 . OOOOO I 

7 . 00000 f 

8.00000 I 

9 . OOOOO ~T 

10.00000 I 
ll.COOOO 1 

12.00000 1 

l3.0t''O'» 1 

M. OOOOO 1 
t:'..00000 I 


• 4 « 

#4 f 



4 4 

4 4 

4 4 

4 4 

4 9 

4 


ITfT I Ifn; IKUllOUl MITH SUITLllO SIOHAS. ITf*t 2 fOR f-tOf 01 ICI U VT. ftR RNS |TIX 3 TO SfOP 

/) /V i ' Figure IV, 19. A Sample Output of APPDEL. 



raiHTlNG C«MM (NINUTES) 










SOIARSIH Subroutine HPTS 


Purpose ; This subroutine computes the power transfer efficiency of 
the spacetenna by calculating the actual power received by the 
rectenna when the beam suffers deflection due to the phase errors 
In Che phasing system and dividing this power by the power received 
by the rectenna when the main beam is Ideally pointed. 



M(I) - Nianber of amplifiers In the 1— level subarray, I « 1,...,10. 

N(I) - Number of radlaging slots in the 1— level, I ■ 1,.,.,10, 

NA(I) - Number of amplifiers In the i-~ level, I ■ 1,.,.,10. 

R{I) - Radius of 1~ power density ring, I • 1,...,10. 

SI6MAB<D-The total standard deviation in the 1~ power density 

level radiating elements, I ■ 1,...,10. 

SI6PHI(I),SI6PSI(I) * The standard deviation of location errors of 

radiating and the receiving elements for the subarrays 
th 

In the 1 — power density level, I « 1,...,10. 

SI6MAI(0) - The standard deviation of the current amplitude feeding 

the radiating elements In the 1— power level, I » 1,...,10. 
CURNT(I) - The current output of the power amplifier in the i~ 
level, I • 1,...,10. 

XhCAN(I) ,YMEAN(I) - The x and y components of mean subarray tilt 

In the 1— power level, I • 1 10. 

S6MA(I) - The standard deviation of the x and y components of the 

subarray tilt in the 1— power density level, 1 ■ 1 10. 

THETAR,PHIR - The (6,;^) direction of the Incoming pilot signal. 
NX,NY,XA,XB,YA,YB,IA,IB - The parameter values of the double 


integral subroutine. 






